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PART I - SMALL DISTURBANCES 


SUMMARY 


An analysis is made of rotating stall in compressors of many stages, 
finding conditions under which a flow distortion can occur which is steady 
in a traveling reference frame, even though upstream total and downstream 
static pressure are constant. In the compressor, a pressure-rise 
hysteresis is assumed. Flow in entrance and exit ducts yield additional 
lags. These lags balance to give a formula for stall propagation speed. 
For small disturbances, it is required that the compressor characteristic 
be flat in the neighborhood of average flow coefficient. 

Restalts are compared with the experiments of Day and Cumpsty. If 
a compressor lag of about twice that due only to fluid inertia is used, 
predicted propagation speeds agree almost exactly with experimental 
values, taking into account changes of number of stages, stagger angle, 
row spacing, and number of stall zones. The agreement obtained gives 
encouragement for the extension of the theory to account for large 
amplitudes. 


INTRODUCTION 


In recent years experimental studies and reviews by Greitzer, Day 
and Cumpsty have clarified and redefined the problem of rotating 

stall of axial flow compressors, giving a new emphasis to the compressor 
of many stages. At low flows, virtually all compressors exhibit zones 
of stalled flow which propagate tangentially, with blade passages 
cyclically stalling and recovering in the manner first described by 
Emmons (^) . 

However, when the compressor has many stages (about three seems 
to suffice) , the stall patterns seem to be simpler and presumably more 
predictable than when the compressor consists of a single row or stage. 
Early reviews (^) convey this impression, but it remained for Day, 
Greitzer and Cumpsty (^) to describe what one might call the kinematics 
of rotating stall. They demonstrated, for example, that the stall 
pattern in "deep stall" tends to run straight through the machine, 
with an especially abrupt transition at the trailing edge of the stall 
zone. 


Theoretical studies in the past have attempted to understand 

rotating stall for a single stage under the natural impression that if 
one could properly relate the Emmons type of stalling-unstalling 
episodes to the individual cascade characteristics, then one could 
proceed to synthesize a theory for multistage machines. This approach, 
though logical, has not been fruitful. In effect, rotating stall is 
too universal a phenomenon and depends on too many variables; single- 
stage analysis does not yield critical conditions meaningful for whole- 
compressor performance, especially in deep stall. 



In the present three-part study, we explore the idea that rotating 
stall theory should begin with the compressor of many stages, operating 
in stall. Doing this, we will not emphasize the problem of "stall 
inception," but, rather, try to incorporate in the theory the kinematic 
features of deep stall described by Day, Cumpsty and Greitzer (_3) . 
Greitzer and Cumpsty (^) have presented a semi-empirical analysis of 
deep-stall propagation speed which is in this same spirit. 

In Part I, a single dynamic hysteresis is assumed for the blade- 
passage force function and conditions are derived for the occurrence 
of a weak stall pattern which is steady in some rotating frame of 
reference. From this analysis we find a formula for stall propagation 
speed which agrees well with experimental data, if the single lag 
parameter is properly chosen. 

In Part II, the analysis is extended to describe the stall patterns 
of arbitrary velocity amplitude which can occur when the quasi-steady 
compressor characteristic is either flat (deep stall) or else curved 
as it would be at stall inception. Wave forms for velocities and 
pressure are found, as well as propagation speeds and the distortion 
of the compressor characteristic due to rotating stall. 

In Part III, a compressor characteristic is adopted which has the 
realistic feature of an unstalled region with peak pressure rise, a 
deep-stall valley, and a steep resistance to reversed flow. Rotating 
stall then is found to have the nature of a limit cycle, familiar in 
the classical van der Pol problem of nonlinear mechanics ( 10 ) . Results 
include a prediction of stall recovery as flow coefficient is increased. 

The present Part I, though limited in results to weak disturbances, 
will include certain developments in forms suitable for more general 
application in the subsequent parts. 

Notation is defined, for all three parts, in Appendix F. 

THE COMPRESSOR 


The compressor is assianed to have a constant-speed characteristic 
as sketched in the solid line of Fig. 1.1. We will denote the upstream 
total (p-p) to downstream static (Ps) pressure-rise coefficient by 
in the absence of rotating stall: 

p -p 

= Hj (4>). ( 1 . 1 ) 

pU ^ 

The axial flow coefficient <{) has the average value 
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( 1 . 2 ) 



The characteristic ipc will be considered to be single- valued, as 
shown. Throttle settings, indicated by dashed parabolic lines, will 
intersect the characteristic and establish operating points (A-E on the 
sketch) about which rotating stall might occur. Points B and D in 
the sketch would indicate limits of a quasi-steady hysteresis loop 
(B-C-D-E) due (we assume) purely to the way the throttle parabolas 
intersect the characteristic. In parts II and III, we shall permit the 
net pressure rise (¥) to depend on the existence of rotating stall. 

For example, the operating point in rotating stall for the $ indicated 
at point A might actually lie above or below the characteristic curve. 

The compressor has N stages, with IGV and OGV rows providing zero 
incidence at design. Blading and design vector diagrams are symmetrical 
(50% reaction) as sketched in Fig. 1.2, Radial effects will be neglected, 
and flows will be considered incompressible. 

FLOW DISTURBANCE IN THE COMPRESSOR 


First, a disturbance will be postulated to exist in the compressor, 
steady in a circumferentially moving frame. This disturbance implies 
disturbances in the flow fields entering and leaving the compressor - 
The plan of this study is to find the requirement that such a set of 
disturbances exists, even though pressures far upstream and downstream 
remain constant; in effect, to find eigensolutions and eigenvalues for 
the whole system. We begin by describing the disturbance of a single 
blade passage. 

Disturbance in a Blade Passage 

^Vhen rotating stall occurs, flow velocities will change locally 
and in time. We will assume that when axial flow is slowly reduced, 
keeping U constant, flow fields of rotor and stators remain symmetrical, 
though pressure rise and deviation should change. In other words, 
the quasi-steady velocity triangles remain symmetrical, and correspond- 
ing pressure changes are equal. Pressure rise in each row is assumed 
to depart slightly from the quasi-steady, however. 

The static pressure rise through each row will be assumed to have 
the form 



= F(<})) - T(c}))i 


whether rotor or stator. 


(1.3) 
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The quasi-steady pressure-rise coefficient F might have the form 
sketched in Fig. 1.3. Assigning the same F to each row, and properly 
taking account of dynamic pressure, F((f)) will yield the of Fig. 1.1. 
The occurrence of a maximum of F, and the drop of F as ({) is reduced, 
will signify stall of a blade passage. 

The second term describes hysteresis due either to flow accelera- 
tion, flow separation or viscous processes in the blade passage. It 
is written in the sense that if <p ±s decreasing, then there is an 
extra force which is positive. At any flow coefficient, therefore, 
we assume a counter-clockwise dynamic hysteresis loop, as sketched in 
Fig. 1.4. The parameter T will generally be a function of (J), of 
course. However, throughout this study, we will take T to be constant, 
and the same for rotors and stators. 

Because Eq. (1.3) will be considered to apply equally to the rotor 
and stator rows, in effect we disregard centrifugal and secondary-flow 
phenomena. In fact, for S 3 mimetric blading, T would be the same for 
rotors and stators. For guide vanes, we will adopt the same equation, 
using subscripts to distinguish them from the stage elements. 

Stall Pattern of the Compressor 


If, by symmetry, rotor and stator blade passages show no deviation 
under disturbance, and no flow development takes place between rows 
(as if there is zero inter-row spacing) , then continuity will require 
that all instantaneously-connected blade passages have the same axial 
velocity. Therefore, any disturbance of axial velocity would be purely 
circumferential, however many stages there are, and whatever the speed 
of propagation of the disturbance pattern. Further study might show 
that dynamic effects will require modification of this symmetry argument. 

Under present assumptions, a stall pattern, shown shaded on Fig. 

1.4, must proceed straight through the machine. This is consistent 
with observations of deep stall. Obviously, if reaction is not 50% 
or blades are of unequal axial length, if compressibility or centrifugal 
effects were significant, or if circumferential flow between rows is 
allowed, then the straight-through pattern would not be correct. 

Next we assume that a disturbance pattern is steady in some rotating 
frame of reference. If the compressor is "unwrapped" as in Fig. 1.4, 
the disturbance is moving to the right with speed fU relative to stators. 
The "stall propagation speed" f is unknown; it is to be found as an 
eigenvalue of the disturbance, so to speak. Relative to the disturbance, 
stators are moving to the left with speed fU and rotors are therefore 
moving to the right with speed (l-f)U. The axial flow coefficient, in 
this frame, takes the form 

(}) = 4) + g(0) (1.4) 
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where g represents a steady disturbance. Again referring to Fig. 1.4, 
distance x is measured to the right, so that 

e = 2x/D (1.5) 


where D is the mean diameter of the compressor. If the disturbance 
of (p is purely circumferential, then Eq. (1.3) applies equally to each 
blade row, and pressure-rise may be accumulated through the compressor 
purely as a function of 9, or equivalently, x. 


9 

The time derivative (p must be expressed in terms of the steady 
disturbance described in Eq. (1.4) and Fig. 1.4; it signifies the 
rate of increase of (p experienced by a blade passage as it traverses 
the postulated stall disturbance. Thus, from Eqs. (1.4) and (1.5), 

} = g’(B)e = g'(6) . (1.6) 


Fig. 1.4 sketches the situation: The solid line represents the steady 

but disturbed flow, (J)(9). Stator blades moving the the left at speed 
fU experience an increase of <{> depending on the local gradient of <p 
(that, is, (})’ (0)) and the relative speed x = -fU. Similarly, a rotor 
moves in the same gradient with speed x = (l-f)U. For stators and 
rotors separately, 


stator 

(p 

rotor 


= -g'(0) £ 

' g'(e) ^ (i-f). 


(1.7a) 

(1.7b) 


Since it is intended to accumulate forces through the compressor, we 
will find use for the combined pressure-rise expression for a single 
stage : 

= 2F((f) - Tg'(6) ^ (f - f> . (1.8) 

stage 



Guide Vanes 


The OGV is assumed to have no deviation, and will therefore deliver 
flow in a purely axial direction to the diffuser. A hysteresis term 
(Eq. (1.3)) applies, with a value of Tqv presumably different from the 
T applicable to stators and rotors of the stages. 

It will be assumed that IGV’s are purely axial at entrance. When 
rotating stall occurs, there will be (as we shall find) a transverse 
velocity disturbance at the compressor (IGV) entrance. In this paper, 
we will consider this velocity disturbance small compared to V. Regard- 
ing the IGV as a nozzle, its entering velocity will therefore depart 
from the local value of v only to second order in small quantities. 
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Accordingly, although there must be a rather abrupt pressure rise 
crossing the entrance plane of the IGV’s, it too will be of second 
order in flow deflection. By our previous stage symmetry argument, we 
should expect the IGV discharge angle to be the same as the OGV entrance 
angle. Thus, recalling Eq. (1.3), there should be a negative value of 
F for the IGV nozzle, just equal to the positive value for the OGV 
passage acting as a diffuser, even though varies. Again a hysteresis 
term should apply, and, for OGV and IGV together, one has, correspond- 
ing to Eq. (1.8), 

i ^ f 

\ 2 / guide 

vanes 

2 

where K^-h represents the coefficient of abrupt pressure rise at the 
IGV entrance. 


UPSTREAM DISTURBANCE 


Far upstream, we imagine a reservoir of total pressure P' 3 ;, constant 
whether or not rotating stall is progress. Between that reservoir 
and the compressor face is a flow passage of some possibly complex 
shape ("Gooseneck"). In that region, a flow disturbance will develop, 
if rotating stall exists. In this study, we shall assume the entrance 
region to be a straight, uniform channel. In that case, if the flow 
upstream has uniform velocity and pressure, 

P<» = Pt - i 

and V will be the average velocity of any disturbance that arises 
downstream. 


Velocity Field 

In the channel, disturbances will be considered irrotational as 
well as incompressible. Thus, Laplace’s equation applies. For con- 
venience, we adopt the coordinate system of Fig. 1.4, so that the 
velocity components far upstream (y = - °°) are 


CO 00 y _ 

~u” IT "" u " 


At the compressor face, Eq. (1.4) applies: 

^ = (j) = $+ g(x) . 


( 1 . 11 ) 


( 1 . 12 ) 
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The resulting steady streamline pattern in the entrance duct region will 
be as sketched in Fig. 1.5. It will be necessary to know the pressure 
disturbance at the compressor face, and for this, the value of u there 
is needed. We may write a disturbance velocity potential 


with $ vanishing when y = - and at y = 0, 


(1.13) 




(1.14) 


where g is arbitrary but periodic, with wave length irD/n. The wave 
number n is introduced to account for the possibility of multiple-cell 
stall. The required u-velocity at the compressor face will henceforth 
be designated as h(x) ; 




(1.15) 


In general, h is the Hilbert transform of g: 


h(e)=-i 


(1.16) 


If g is periodic, it may be shown that Eq. (1.16) yields 


2tr 


h = - ^ / g' (C) In I sin 


0 -^ 


dC. 


(1.17) 


If g(0) has a Fourier series. 


g(0) = I (n sin n0 + b cos n0) 
" n n 


(1.18a) 


then 


2y 


$(y,0) = I e ^ (a^ sin n0 + b^ cos n0) 


(1.18b) 


which satisfies Eq. (1.14) and vanishes at y = - The corresponding 
u and v velocities both vanish as they should at y = - «», and 


h(0) = y (a cos n0 - b sin n0) 
" n n 


(1.18c) 
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For example, if g is a harmonic wave, 

g = sin n0; h = cos n0, (1.19) 

£ being an amplitude parameter. 

Pressure Field 

In the disturbance coordinates, the steady Bernoulli equation 
applies; 


p^ + I pU^(fV$^) = Pg + J pU^[(-f+h) V(<IH-g)^] 


( 1 . 20 ) 


where is the pressure at the compressor face, just ahead of the IGV. 
Recalling Eq, (1.10), one finds 

12 12 12 

= -(^ r + $g + y - fh + ^ h ). (1.21) 

pU 

If g is a weak harmonic wave, Eq . (1.19), with e assumed small 
compared to 1, would provide 


Po'Pt 


PU 


12 2 
Y $ - e($ sin n0 


- f$ cos n0) 


to first order in c. 


( 1 . 22 ) 


If the channel is 
not irrotational, Eqs. 


not straight, or if the 
(1.21) and (1.22) would 

DOWNSTREAM DISTURBANCE 


flow in the channel is 
require modification. 


Far downstream, we imagine a flow of uniform static pressure (Pg) » 
but not necessarily constant velocity; a shear flow might be the final 
result. The flow is expected to emerge from the OGV in an axial 
direction (viewed in stator coordinates) but with varying pressure (pe) 
and perhaps vorticity. The actual exit-flow device, or diffuser, is 
very important in this regard, and we will limit this discussion to 
two types, a straight channel, and a sudden expansion. These two 
extreme types perhaps bracket the possibilities. 

Straight Channel 

If there is a straight, constant-area channel downstream of the 
compressor exit (Fig. 1.6a) then the Euler equations must be solved in 
the region sketched in Fig. 1.5. Defining a stream function for that 
region. 
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(1.23) 


— = l|;; — =(J) = -i|; 

U * U ^ 


and expressing pressure in the dimensionless form 


= P (1.24) 

PU 

where pg is the (constant) downstream static pressure, one may write 
the Euler equations as 


-P = ~ Ip ip 
X y xy 


]p ip : 

X yy 


-P = iL - ip ip 
y y XX X xy 


(1.25) 


Boundary conditions at the rear face of the compressor (subscript e) 
at y = 0 are 


and far downstream (subscript s) at y = <», the pressure must become 
uniform, but vortical velocity disturbances may survive: 

Pg = 0. (1.27) 


In formulating a solution of the system of Eqs. (1.25-27), we 
should recall that only P is needed, in its relation to the disturbance 
function g(x) . Eq . (1.25) can be converted to 



-Ip Ip + (Tp ) 
XX yy xy 


and Eqs. (1.25) and (1.26) yield 


(1.28) 




(1.29) 


Of course, one or the other of Eq. (1.25) must be included with Eq. 
(1.28). However, if the disturbance is weak, then the right side of 
Eq. (1.28) is of second order, and Eq. (1.28) is simply Laplace's equa- 
tion 
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V^P = 0. 


(1.30) 


The desired solution for P can be found directly, using Eqs. (1.27) and 
(1.29) as boundary conditions. 

In fact, the function -P/f would correspond to the x-derivative 
of a velocity potential, vanishing at y = °o, of which g(x) is the 
y-derivative at y - 0. But that is precisely the relationship between 
g and h defined in Eqs. (1.14) and (1.15). Taking account of the 
difference of direction of y between the two cases, one may immediately 
conclude that 


P = fh; P 
e e 


P “P 
s e 



(1.31) 


This is not a general result, because the right side of Eq. (1.28) was 
set equal to zero. It is, however, correct to first order in small 
quantities if, for example, Eq. (1.19) applies, and may be a useful 
relation even if disturbances are not weak. 


If disturbances are indeed weak, Eqs. (1.25) and (1.31) may be used 
to find the velocity field in the straight channel. The result for Tp is 


if = I- {_^(e + 2 1 + 


e ^ 2 [f sin n(0 + 2 -I 

$ +f 

+ cos n(0 + 2 -| ^) 
f D 

+ e (-$ sin n0 + f cos n0)J}. 


(1.32) 


If disturbances are of finite strength, and if g(0) can be expressed as 
a Taylor series in e, then the full Eq. (1.28) could be made the basis 
of a method of successive approximations. 

Sudden Expansion 

If there is a sudden large transverse (radial) expansion just at 
the OGV exit, as sketched in Fig. 1.6b, then the exit flow would be 
like a jet emerging from a long slot. While the jet velocity would 
vary circumferentially on a scale comparable to wheel diameter, the 
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pressure matching of the jet to the receiver would occur on the scale 
of blade height. Since the receiver pressure could be considered 
constant, the jet exit pressure could also be considered constant. In 
this extreme of an orifice-like diffuser. 


= P. 


( 1 . 33 ) 


This is the simplest possible assumption about pressure change at the 
compressor exit. Of course, it envisages an extremely short diffuser, 
which would generally not be acceptable because of its low efficiency. 

Apparently, the general problem of finding the pressure change at 
the compressor exit is quite difficult, both to define and to solve. 

The specifics of diffuser configuration would greatly affect the result. 
Here, we have found expressions for two extreme cases which may possibly 
bracket t 3 rpical practical situations. We may, for convenicence, write 
them in the form 


P 

e 


P -P 
s e 



(m-1) fh 


( 1 . 34 ) 


where m = 1 would signify a sudden area change, and m = 2 would signify 
a straight channel. We will regard m, therefore, as a somewhat arbitrary 
parameter of our problem. 

REQUIREMENTS FOR ROTATING STALL 


The various pressure relations may now be combined to express a 
statement that pressures remain constant in upstream and downstream 
reservoirs, even though an unsteady disturbance occurs in the compressor 
which connects them. We have already specified that the disturbance 
is steady when viewed in a frame rotating with velocity fU. 


Pressure Match 


We require that Pg“P.p be constant, or 


Po'Pt . Pe"Po 

r\ "r r\ 


pu 

entrance 

duct 


pu 


s e 


pu 


Ps"Pt 


pu 


E m . 


compressor exit duct reservoirs 


( 1 . 35 ) 
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The first term is given by Eq. (1.21), the third by Eq. (1.34), 
and the second by application of Eq . (1.9) for guide vanes and Eq. 
(1.8) for each of N stages. The result may be written 


. 2UNT 
^ D 





f]g*(e) - mfh h^ 


- I Kjh^ + m - [NF((|)) - I (J)^] = 0. 


(1.36) 


It is obvious that the last quantity in brackets is the compressor 
characteristic function that would apply in the absence of rotating 
stall; that is, the curve i|^Q((f)) shown on Fig. 1.1. 

ip^(4>) = NF((j)) - j (f)^. (1.37) 

It should be emphasized that while 4^c(4*) is a given feature of the 
compressor, it is not the same as 'F; the latter is the pressure rise 
actually produced, with rotating stall (if it occurs) in progress. 

While ipc is, in principle, given, it is not the characteristic usually 
measured in tests of compressors subject to rotating stall. One might 
rather say that ¥ is the usually-measured quantity; here, it is part 
of the answer. 

The problem now is to find what values of parameters such as f, m, 
'{', or features of the function i|^^(4)), permit nonzero solutions of the 
pair of equations (1.36) and (1.37). 

Weak Disturbances 


In this paper, we will restrict disturbances to be weak. Thus, 
we neglect the terms (l/2)h^ and (l/2)Kih2 appearing in Eq. (1.36). 

Also, since (f) will remain close to $ during a weak disturbance, we 
keep only the first two terms of the Taylor series 

+ ... (1.38) 

For convenience we define 
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f) - 


f 


(1.39) 


, _ 2UNT .1 
^ “ D ^2 


U(t +T ) 
^ IV OV'^ 


so that Eq. (1.36) becomes 


Ag’(0) - mfh - + 'P - ii>^($) = 0. 


(1.40) 


By definition, the average value of g over a cycle must vanish. 
Also, if the flow far upstream is purely axial, no average swirl can 
arise between there and the compressor entrance if the flow is irrota- 
tional. Therefore, the average value of h must also vanish. These 
conditions would apply even for strong waves, but in the present case, 
integrating Eq. (1.40) over a cycle shows that 


4' - ij; ($) = 0. 

c 


(1.41) 


In other words, for weak waves, there is no vertical displacement of 
the operating point from the rotating-stall-free characteristic; it 
is as indicated for point A of Fig. 1.1. 

Presumably, g and h have Fourier series with coefficients related 
by Eqs. (1.18a,c). Because Eq. (1.40) is linear, each coefficient is 
independently fixed. Using only the a^ terms, we find by substitution 


(Xn-mf)a cos n0 - iIj’($) a sin n0 = 0. 
n ^ n 


(1.42) 


The conditions for a weak periodic wave of arbitrary amplitude are 
therefore that Xn = mf, that ijJc('^>) = 0, and that = i|;c($) . The first 
of these may be solved for propagation-speed coefficient f: 


f 


1/2 

1 4 . D . 2 
^ ™ 2UnNT Nt 


(1.43) 


This is the chief result of our analysis, and will be discussed fully 
in the next section. The second requirement 
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(1.44) 


= 0 


means that the local curve of ipj, must be flat. This could be true in 
the deep-stall region (points like A or C in Fig. 1.1), and is true 
at the maximum of the characteristic (point M) . In the latter case, 

Eq. (1.44) would be a condition for stall inception, agreeing with the 
classical stability requirement. On the descending branch of the speed 
line, Eq. (1.44) would not be satisfied; we interpret that to mean not 
that rotating stall cannot occur, but rather that it cannot be weak. 

In fact, ip^ can be interpreted as a negative damping factor; since Eqs. 
(1.18) indicate that g = -h’, Eq. (1.40) is that of a harmonic 
oscillator. 


Xh" - i|^'h' + mfh = 0 
c 

with homogeneous boundary conditions, a wave must decay if ipc is nega- 
tive, or be amplified if is positive. In neither case could a 
permanent oscillation exist. 

Eqs. (1.41), (1.43), and (1.A4) therefore give the conditions for 
a steady rotating stall of weak amplitude to exist, when reservoir 
pressures are both constant. 

STALL PROPAGATION SPEED 


We next discuss the implications of Eq. (1.41) and compare its 
predictions with experiment (_2) . Eq. (1.43) is essentially a balance 
of disturbance phase lags inside and outside the compressor. 

The External Lag Parameter 


The parameter m appears in Eq. (1.43). It was defined in Eq. 
(1.34) as a measure of the strength of the lagging pressure term fh 
arising in the exit duct. The same lag term arises from the entrance 
duct (the term fh in Eq. (1.21)) assuming ideal flow in a straight 
channel. We should admit the possibility that in the entrance flow, 
too, configuration or loss effects might affect lag. Thus, where the 
two lags are added to give mfh in Eq. (1.36), we should imagine m to be 
a net lag parameter depending on the combined contributions of both 
entrance and exit flows. We do nevertheless expect that it depends 
chiefly on the exit arrangement. 
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The Compressor Lag Parameter 


The lag T, defined in Eq. (1.3), depends on the aerodynamics of a 
blade passage, sketched in Fig. 1.2b. Whatever its physical nature, it 
should scale with axial passage length L and axial velocity V. We are 
then led to define a dimensionless lag 


T* 


_ $UT 
~ L 


(1.45) 


in terms of which Eq. (1.43) becomes 


f 


1/2 

T* 

m$ D _1 _1 __v 

n 2NL T* N T* 


(1,46) 


The lags for guide vanes are assumed equal. 


It is perhaps of interest to note the value of T* that would 
represent "apparent mass" or inertia only. Referring to Fig. 1.2b, 
we see that a mass pLa/cos y of fluid is accelerating in a blade 
passage at a rate U(|i/cos y requiring the force provided by an additional 
pressure rise (Ap)^: 



2 


2 

cos y 



(1.47) 


Comparing Eq. (1.3) we see that, due to inertia in blade passages only, 
T* = 2$/cos^y. The stagger angle y is roughly tan“l(l/2 4><^) where 
is design flow coefficient. We write 


2 $ 

T* = - ^ k (1.48) 

cos y 

where a factor k has been included because, often, the interrow gap is 
quite large. If the gap were a multiple k-1 of L, then T* would be 
proportional to k. 

Other effects having to do with separation or secondary flows will 
also influence T*. If these could be represented by suitable choice of 
k, it would be appropriate to write Eq. (1.46) in the following form: 
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1/2 


(1.49) 


f 


m , 2.1 

1 + ^ (cos Y) 


T* 1 
V 1 


2n 2NL T* N 


We will use k as a general "uncertainty factor" comparing actual T to 
blade-passage inertia lag. It is interesting that the external and 
internal lag parameters enter Eq. (1.49) only in their ratio, m/k. 

Qualitative Results for Propagation Speed 

Eq. (1.46) shows certain general trends: 

1. Propagation speed is less than 50 percent of U, unless T is 
negative. Actually, negative lag has been argued to occur at stall of 
isolated airfoils (11) , and f > 1/2 is generally associated with single- 
stage experiments (^»^) • Perhaps these results are consistent. 

2. Increasing number of stages (N) increases f toward 1/2. This 
result has been found by Cumpsty and Greitzer (^) . It is important 
that N multplies T* in Eq. (1.46). If T is very small, being propor- 
tional to blade chord, it might at first be thought to make a negligible 
contribution to pressure rise (Eq. (1.3)), and therefore not be proper 
to include in a disturbance analysis. This argument might be valid for 
single-row problems. However, a small lag per row will accumulate over 
many rows in a multistage compressor to make a large overall lag; it 

is the latter which governs rotating stall. The same idea is expressed 
by noting the group D/2NL in Eq. (1.49); it describes the shape factor 
of the whole compressor when viewed from the side. Typically, it would 
be of order 1, while for a single row, D/L would be much larger — 
typically about 15. 

3. Multicell stalls have higher propagation speeds. The mode num- 
ber n represents number of stall cells. 

4. Increased lag (larger k) tends to increase propagation speed. 

5. Especially if design flow coefficient is small, smaller 
(or larger stagger angle) will increase f. 

6. If the net external phase-shift parameter ra is made smaller, 
f is increased. For a typical engine, m is probably near 2. If the 
diffuser is a sudden area change at the rear face of the compressor, 

m = 1. If circumferential velocity were suppressed in the inlet flow, 
m might also decrease by 1. A value of m = 0 may be theoretically 
possible, if circumferential velocity is suppressed in both entrance 
and exit flows . 


7. Guide vanes have little effect on f in multistage compressors. 
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Comparison with Day and Cumpsty’s Data 


The foregoing results may be compared with the propagation speed 
data of Day and Cumpsty (^) . Table 1.1, in the first 5 columns, 
collects some of their results for various designs, number of stages, 
and number of stall cells. Of single-cell cases (n = 1), only the 
full-span results are included, because the present paper assumes full- 
span stall, in effect. Multiple-cell cases are included, but it should 
be noted that all of them are part- span. Guide vanes have a smaller 
average incidence than core blades, by about a factor of 2. Therefore 
relative velocities are lower. As a rough estimate, we will assume 
to be the average of Eq. (1.48) and 2$k. In that case, 

T* , ^ 

For the experiments (^) the blade chord is 17.8 mm, and the mean 
diameter is about 320 mm. The spacing between rows is about one chord 
width. In that case, if inertia were the only source of lag, we would 
expect k = 2. There appears to be a fairly long straight section 
before the diffuser, and so one would expect m = 2. Eq. (1.49) shows 
that k and m enter in the ratio k/m. The last two columns of numbers 
are the theoretical predictions of Eq. (1.49) for two values of k/m, the 
first being for the expected k/m =l(k=2,m=2), and the second 
for k/m = 2. The latter would apply if k were 4 instead of 2 owing 
to separation hysteresis, or perhaps if the exit acted like an orifice 
so that m = 1 instead of 2. The former would be a more likely combina- 
tion (m = 2, k = 4) . 

Comparing the theoretical and experimental columns, it is clear 
that Eq. (1.49) gives all the proper trends of f with Y» N, and n. 

We should keep in mind that the theory does not contemplate the part- 
span stall patterns found experimentally for multiple zones, so agree- 
ment for n different from 1 is perhaps fortuitous. The theory 
probably does not apply well to single-stage situations. Therefore, 
the most appropriate comparisons are for n = 1, N = 2,3,4. These 
cases are enclosed in boxes in the Table, for emphasis. 

Table 1.1 clearly indicates that k/m = 2 is the best choice for 
agreement, and the agreement could hardly be closer, for intermediate 
and high-$d builds. For low full agreement would require a some- 
what higher value of k/m, about 3 instead of 2. Perhaps the higher 
stagger angle causes more hysteresis due to separation. 

In the low-'fjj experiments, a reduced spacing gave a substantailly 
lower speed. In that case, we would expect to reduce the factor k to 
account for the smaller mass of fluid undergoing oscillation. Because 
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spacing was originally one chord length, it seems reasonable to reduce 
k by a factor of 2. Indeed, Table 1.1 shows that k/m = 1 gives a 
prediction having the same relation to experiment for reduced spacing 
as k/m = 2 had for the wider spacing. Also, it was mentioned that 
k/m = 3 gives exact agreement with experiment for the wider spacing; 
k/m = 1.5 gives exact agreement for reduced spacing. Thus, we may add 
spacing to the list of compressor parameters whose influence on propaga- 
tion speed is well predicted by Eq. (1.49). 

Finally, we should note that Cumpsty and Day (^) have derived an 
equation for propagation speed based on a consideration of pressure 
jumps across the boundaries of deep-stall zones, together with certain 
observations of compressor test results. Their equation also predicts 
the correct trend with number of stages. 

CONCLUDING REMARKS 


The present analysis shows that rotating stall can be considered 
an eigenfunction of a compressor- inlet-diffuser system if a single lag 
process occurs in the generation of force by a blade passage, and, at 
least for weak amplitude, if there is a flat total-to-static pressure 
rise characteristic in the neighborhood of the average flow coefficient 
set by the throttle. 

The predicted propagation speed depends on various features of the 
compressor, and follows very closely the trends of experimental observa- 
tions as those features are changed. The lag processes originating 
in the compressor (k) and those originating outside the compressor in 
the entrance duct and the diffuser (m) require more study from a funda- 
mental point of view. 

The present analysis, being based on small perturbations, has 
limited value. First, actual perturbations in deep stall are not 
small, nor are amplitudes arbitrary. Therefore, there is little 
guidance from the present analysis as to how compressor parameter 
changes might help to control or eliminate the phenomenon. However, the 
results of the present analysis perhaps do give encouragement for the 
study of finite amplitude oscillations using similar physical assumptions. 

Parts II and III of this study report such extensions. In Part 

II, we keep the restriction to a flat characteristic, but allow the 
amplitude to be finite. Results will include a measure of how the 
mean performance is affected by the existence of rotating stall, as 
well as certain interesting velocity and pressure profiles. In Bart 

III, we attempt to relate rotating stall to the complete characteristic 
curve. Results include an indication of when the "blockage" model (^) 
applies, and also predict that rotating stall becomes impossible beyond 
some value of flow coefficient, probably indicating the onset of stall 
recovery . 
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TABLE 1.1. PROPAGATION SPEED COEFFICIENT 


Y (degr) N 

n 

50 1 

1 


6 


12 

50 3 

1 


4 


12 

50 (reduced ^ 

1 

spacing) 

± 

5 

35 1 

1 

2 

1 

3 

1 

4 

1 

20 1 

1 

2 

1 

3 

1 

4 

1 

Chord = 17.8 mm 


Ref 2 


f (percent) 

Eq. (1.49) 
k/m = 1 k/i 


Mean diameter: 320 mm 
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Fig. 1.1. Constant-speed compressor characteristic operating 

points ; upstream total to downstream static pressure 
rise versus axial flow coefficient. 



(o)Sketch of developed (b)Sketch of single 
compressor blade passage 


Fig. 1.2. Sketch of compressor configuration. Number of stages is 
arbitrary; two are shown. Radial effects are neglected. 
Entrance and exit planes are denoted by @ and @ . Blading 
is symmetrical. 
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Fig. 1.3. Assumed form of blade-passage pressure-rise function. 



Fig. 1.4. Sketch showing axial stall zone in compressor, in coordinates 
fixed in the stall zone which, in laboratory coordinates, is 
moving to the right at speed fU. 


21 




compressor 




XorS 


Fig. 1.5. Sketches of flows in straight entrance and exit ducts, 

coordinates fixed in disturbance. At the compressor exit, 
u = -fU. 




Fig. 1.6. Sketches of extreme exit-duct configurations. 
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PART II - FINITE DISTURBANCES 


SUMMARY 


A small-disturbance theory of rotating stall in axial compressors 
is extended to finite amplitude, assuming the compressor characteristic 
is a parabola over the range of the disturbance. An exact solution 
is found which requires the operating point to be at the minimum or 
maximum of the parabola. 

If the characteristic is flat in a deep-stall regime, the previous 
harmonic solution applies with neither reverse flow or "unstalling". 

If the characteristic is concave upward in deep stall, the disturbance 
has a skewed shape, steeper at the stall-zone trailing edge as experi- 
ment shows. Propagation speed is only slightly affected by this non- 
linearity. 

Near stall inception, negative curvature in combination with 
multiple stall zones can limit the nonlinear oscillation, in the manner 
of "progressive" stall. If, as seems likely, lag at stall inception 
is negative (opposite to inertia) , propagation speed exceeds 1/2 wheel 
speed, as experiments suggest. 

INTRODUCTION 


In Part I ("Small Disturbances") , it was shown how a combination 
of three lags, pertaining to entrance-flow passage, compressor blade 
passage, and exit-flow passage permit the occurrence of a flow distor- 
tion which is steady in a moving frame of reference. In principle, the 
entrance and exit lags are easily calculated, while the blade-passage 
lag is difficult to define and infer from tests. Comparison with 
measurements of rotating-stall propagation speed made by Day and Cumpsty 
(^) support the idea that the blade-passage effect is largely inertial, 
in deep stall. 

The lag per row was represented as a factor k of the purely inertial 
value for a blade passage. If the inter-row gap is equal to the blade 
depth then k would be about 2, still considering only inertia but 
including the mass of fluid occupying the inter-row space. The combined 
lag effect of entrance and exit flows was represented by the factor m; 
if long straight passages describe both entrance and exit configurations, 
then m = 2. 

For small perturbations, k and m enter as a ratio in the formula 
for propagation speed. If the ratio is taken as 2 (about twice what 
we would expect on purely inertial grounds) , agreement with experiment 
is very complete, suggesting that a single compressor lag parameter, if 
suitably chosen, could lead to a quite complete description of rotating 
stall. 
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Small perturbations restrict wave forms to be harmonic, which is 
quite unlike the wave forms actually observed (^) . Furthermore, the 
amplitude is unspecified in a linear analysis. Thus, although a good 
prediction of propagation speed is encouraging, the important questions 
about rotating stall have to do with disturbances of finite amplitude. 
Specifically, 

1. What is the effect on average pressure-rise performance caused 
by rotating stall? 

2. Are there circumstances which limit the amplitude of oscilla- 
tion? If so, "progressive stall" may be indicated. 

3. Are there ranges of flow coefficient for which rotating stall 
is forbidden? If so, "stall recovery" may be indicated. 

In considering these questions, we should, if possible, allow the 
assumed compressor characteristic to depart from a constant value during 
the oscillation of flow coefficient. In this part, we investigate 
the finite-strength disturbances permitted when the characteristic 
is either flat or has constant curvature. 

SOLUTIONS OF THE PRESSURE EQUATION 


The Disturbance Equations 


We recall the basic pressure-match equation, Eq. (1.36) of Part I, 
and incorporate Eqs. (1.37) and (1.39): 

Xg'(0) - mfh + ~ (l-K^)h^ + ’i' - iP^((j)) =0 (2.1) 


where g and h are the axial and circumferential velocity disturbances, 
related by the Fourier series (Eq. (1.18)), or the equivalent statement 


h + ig = analytic function of 


(e ) . 


( 2 . 2 ) 


We will assume that the rotating-stall-free characteristic function 
remains well represented, over the range of the oscillation, by the 
first three terms of a Taylor series (Eq. (1.38)) about the mean flow 


ip^(<p) = i/;^(#) + i/jV$)g + I + ... 


(2.3) 
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We will find that no solutions are possible unless = 0, even for 

finite disturbances. Such a term would represent a damping effect, 
even for finite disturbances, and must be discarded if only the 
three terms of Eq. (2.3) are considered. The reasons why this is true 
will be made clear later, especially in Part III. 

Eq. (2.1) now becomes 

Xg’(0) - mfh + I (l-K^)h^ - ~ \l)’^($)g^ + 6 = 0 (2.4) 

where 6 is the performance increase (perhaps negative) due to rotating 
stall : 


6 E H' - i^^($) . 


(2.5) 


A typical function positive 6 are sketched in Fig. 2.1. 

The second term of Eq. (2.4) embodies the factor m, introduced 
in Part I to account for lags external to the compressor. This factor 
was found to be 2 for small disturbances of a straight-channel exit. 

In the present part, we are concerned with finite disturbances, and a 
more complete accounting for downstream pressure recovery is no doubt 
needed. The necessary analysis would be very arduous, and not justified 
for a particular highly-idealized exit channel. For a sudden exit 
expansion, we argued that m = 1, and this result is independent of 
amplitude. Therefore, the term mfh will be retained because it is proper 
when m = 1, and may be empirically justified for other values of m. 
Further research on this point is required. 

The disturbances g and h governed by Eq. (2.4) are periodic (but 
not necessarily harmonic) and have vanishing averages over a cycle; in 
the case of g because it is defined as a departure from average axial 
velocity, and in the case of h, because net circulation in the entrance 
flow is assumed to remain zero. Therefore, integrating Eq. (2.4) over 
a cycle shows that, when nonlinear terms involving h^ and g^ are 
included, the constant 6 cannot be zero, as it was in the linear 
analysis of Part I. 

A General Solution 


In Part I, Eqs. (2.2) and (2.4) had a simple harmonic solution 
because the last three terms of Eq. (2.4) were absent. Here, those 
terms are kept, and appropriate nonlinear periodic solutions are sought. 
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It happens that a general solution, believed to be unique, can be 
written down at once: 


h + ig 


inO 


1 + Tie 


in0 


(2.6a) 


or 


, A(cos n0 + Ti) 

h = 

1 + r| + 2rt cos n0 


g = 


A sin n0 


1 + n + 2n cos n0 


(2.6b) 


where A and ri are constants to be determined by satisfying Eq. (2.4). 

In order that h and g have Fourier series in 0 of the form of Eq. (1.18) 
it is necessary that |n| <1. Eq. (2.6a) is clearly a special case of 
Eq. (2.2) and the previous linear solution is recovered by setting T) = 0 

The forms of these velocity traces are shown for various values of 
q between 0 and 1, on Fig. 2-2. The axial-flow disturbance g is anti- 
symmetric and is traveling to the right. Its negative portion repre- 
sents a "stall zone"; if A is positive, its "trailing edge" is steeper. 

In presenting the results for A and q obtained by substitution into 
Eq. (2.4), we note that f is not known, and recall that 


X = 


2UNT 

D 



2Ut 

V 

D 


f 


(2.7) 


if the lags for IGV and OGV are equal. Thus, a free constant will 
remain after q. A, and 6 are given in terms of f. Results are 


+ (i-5n^Ki-Kj,) 
(i+n^)ip;: + (i-3n^) (i-Kj.) 

4mf 

h 2 2 

’^c (l+q^)r + (l-3q^)(l-K^) 

c I 


( 2 . 8 ) 


(2.9) 
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6 


( 2 . 10 ) 


= ^ [4n - + 1 - K )]. 

8tl^ " ^ 


The next step will be to describe rotating-stall features implied by 
Eqs. (2.8-10). First, we should specify the parameters and K^. 

In these equations, Kj is the coefficient of pressure recovery 
in the entrance of the IGV passage, owing to the fact that a transverse 
velocity disturbance makes an incoming stream tube somewhat narrower 
than the IGV entrance. Fig. 2.3 sketches the process. The conserva- 
tive estimate of this effect is gotten by assuming the area increase is 
sudden. The usual pressure-rise formula in such a case is 


Ap 

2 2 

p(V>u^) 



/ 2 2 

but the geometry of Fig. 2.3 shows that A 1 /A 2 = V// V +u . Since 
V = U(j) and u = Uh in stator coordinates, we may conclude that 


l+(h/(j))^ - 1 ] = 7 hV, 
PU 


( 2 . 11 ) 


the last member of the equation being the definition of Kj. The quantity 
1-Kj which appears in Eqs. (2.8-10) is the head loss coefficient due to 
oblique entrance to the IGV. Clearly, it is a function of h itself, and 
therefore cannot properly be used in Eq. (2.4) as a constant. Fortunate- 
ly, the magnitude of 1-Kj is very small, not exceeding 0.1, by Eq. (2.11), 
unless h/<J) exceeds 0.7. Therefore, one is justified in neglecting 
the head-loss effect represented by 1-Kj, and Eqs. (2.8-10) become: 


Xn 

mf 



A = 


4mf Ti(l-n ) 
2 


ilj" 

^c 


1 +n 


6 


4(mf)^ n^(i-n^) 
f 2.2 ■ 

^c (l+n ) 


( 2 . 12 ) 


(2.13) 


(2.14) 
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We bear in mind that X is a given function of f (Eq. (2.7)) and there- 
fore Eq. (2.12) is an equation connecting f and r|. In fact, Eq. (2.12) 
may be replaced by 


f 


1/2 


1 + 



am 1-Ti 


i+n 


2 

2 


where the symbol a is introduced to define the group 


(2.15) 


_ D _ _D_ $ 

“ 2imx “ ML 2 t* 


(2.16) 


and the last expression results from Eq. (1.45). In discussing compres- 
sors of many stages, one might well neglect the second term of the 
denominator of Eq. (2.15), finding 


f = 


1/2 


1 + 


am 1-ri 

n 2 
1+n 


(2.17) 


Pressure Disturbance 


We should record for subsequent use an expression for pressure- 
rise fluctuation, because that is a quantity more easily measured than 
the velocity components represented by g and h. Eqs. (1.8) and (1.9) 
and the present Eqs. (2.4) and (2.5), with Kj = 1, yield the static 
pressure rise across the compressor: 


C 

P 


P “P 19 

- [i + fCl.) ] 
pu 


-mfh + ^g + Y 


4 . 1 1.2 

! + 2 h ■ 


(2.18) 


Included in brackets on the left of this equation are certain quantities 
which will be constant during any rotating-stall process. Quantities 
which will fluctuate are on the right side, and they are the ones of 
most interest. 
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Amplitude of Solution 


We have noted that Eq. (2.6b) describes an oscillation of axial 
flow (g) that is symmetrical: its amplitude can be found using Eqs. 
(2.13) and (2.17) to be 


2m 




n 


, amv , 2,, amv 

c (1 + — ) + n (1 - — ) 
n n 


(2.19) 


with ri being an amplitude parameter as yet undetermined 


As a function of q, the disturbance amplitude gj^ might be limited 
in some circumstances. According to Eq. (2.19), the maximum value that 
can assume is 


_ _ra 1 

®mm ^ 

^ V l-(am/n) 

when 


n 


mm 


/l-(am/n) 
/ l+(am/n) 


( 2 . 20 ) 


( 2 . 21 ) 


In such a case, Eqs. (2.17) and (2.14) give 


f 

mm 


1/2 

1 - (am/n)^ 


6 

mm 


m am/n 

^c [1 - (am/n)^] 


( 2 . 22 ) 


Because rt is restricted to be less than 1, Eq. (2.21) shows that this 
particular limiting solution is possible only if a is negative, and 
that requires T to be negative; in effect, the hysteresis loop of 
Fig. 1.3 would be clockwise in such a case. 

If, on the other hand, a is positive, Eq. (2.19) shows that maximum 
would occur when ri = 1, taking the value g nnn = 

Next, we should apply these results for interesting choices of 
indicated on Fig. 2.1. There are two types of points where the charac- 
teristic (in absence of rotating stall) is horizontal 

yet has nonzero curvature. They are indicated on Fig. 2.1 as points C 
and M. At point C, the curvature i|;c($) is positive, while at M it is 
negative. 
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DEEP STALL 


Flat Characteristic 


At a point such as C of Fig. 2.1, we first consider the possibility 
that ijic = 0 everywhere in the range of oscillation. In that case, Eq. 
(2,19) shows that D must be zero if g^ is to be finite, and Eqs. (2.17) 
and (2.14) give 


f = m 

1 + (am/n) * 


6 = 0 . 


(2.23) 


This is simply the harmonic solution discussed in Part I, with, accord- 
ing to Eq. (2.6b), 


h = g cos n0; g = g sin n0. 
m m 


We conclude that in a region where the characteristic is flat, 
the solution is the linear one with arbitrary amplitude, even though a 
nonlinear solution is allowed. This result perhaps explains the surpris- 
ing outcome of Part I, that a small-perturbation theory yields a good 
representation of stall propagation speed, even though the actual stall 
amplitude is large. For realistic comparison with experiment, the lag 
T, and hence a, is positive, corresponding to values of f < 1/2. 

The amplitude cannot actually be unlimited, of course. We suppose 
that this linear solution will be appropriate for rather small mean 
flow coefficients, so that the oscillation does not intrude into the 
rising part of the characteristic, with "unstalling" becoming a feature 
of the oscillation. In the same way, we would forbid the flow to 
reverse, imagining an abrupt rise of the characteristic below cj) = 0. 

If this is true, the oscillation should have an amplitude limited by 
<l> itself. That is. 


h = $ cos n0; g = 0 sin n0. 


(2.24) 


Although the velocity components are harmonic, the pressure fluctuation 
is not, according to Eq. (2.18): If n = 1, 


C = -y - mf 0 cos 0 + sin 0. 
p 4 


(2.25) 
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If, for example, $ = 0.4, a = 0.4, and m = 2, the velocity and 
pressure-rise fluctuations would be shown in Fig. 2.4. The pressure 
trace is also harmonic. The propagation is toward the right of the 
figure, at a speed (f) of 0.278. In this case, of course, 6=0, 
because ti = 0. The operating point lies at point C, on the charac- 
teristic, as indicated in the sketch. 

Concave Characteristic 


Next we consider the possibility that is slightly positive, 

or concave, as at point C of Fig. 2.1. Again, we expect that a is 
positive (on the ground that deep stall seems always to have f < 1/2). 

In that case, the greatest possible value of r\ would be 1, but then only 
if is larger than m/g^ according to Eq. (2.19). For realistic 
diagram shapes, must be much less than m/<J> , and we therefore conclude 
that n is quite small in these cases. 

From Eq. (2.19), setting gm = ^ and making r) small gives 


n 


(1 + ^) . 

2m n 


(2.26) 


Substitution of Eq. (2.26) into previous equations will give correc- 
tions of order to the flat-characteristic solutions. While the 
predicted traces of g and h, as well as 6, will change to first order 
in the value of f will be affected only to second order. 

If > Oj then A is positive and the velocity traces are as in 
Fig. 2.2, with an axial flow disturbance (g) which is steepest at the 
stall-zone "trailing edge", as experiments have shown (^) . By way of 
example, we modify the example displayed in Fig. 2.4, by making ijj’^(^') = 
1.65; the changed characteristic is sketched in Fig. 2.5. Eq. (2.18) 
gives the corresponding pressure trace Cp. Pressure rise as well as 
axial velocity demonstrate the steep "trailing edge" found experimentally. 

INCIPIENT STALL 


Convex Characteristic 


At a point such as M, the value of is negative. At a given 
speed, as flow coefficient $ is reduced, rotating stall can first 
appear at point M. The characteristic of Fig. 2.1 is drawn to show 
some degree of overturn to the left of M, terminated in an abrupt drop 
at point B. At that point, we imagine a certain positive slope S to 
occur. The final slope S depends on the compressor; lightly-loaded 
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compressors tend to permit a larger S before a steep break occurs. 

If a parabola is used to represent the curve in the stall-inception 
zone, centered on point M, one could relate S to and the distance 
between points B and M on the diagram: 


A$ = S/ip”. 


(2.27) 


As we have seen, a stall oscillation about point M would have a 
maximum permissible amplitude g^ depending on the sign of t (or a) : 


m/ij;;; 


ram 


, a < 0 


/l - 


(am/ n) 


( 2 . 28 ) 


= "/^c 


, a > 0 


If this maximum amplitude exceeds the A$ between points B and M, then 
the oscillation could easily extend into the deep stall zone, and a 
different type of larger-amplitude oscillation would result, of the sort 
to be discussed in Part III of this report. On the other hand, if gnj^j < 
A<l>, the oscillation could be limited to the vicinity of point M, 
governed purely by the local convex parabolic segment. This latter 
situation would presumably relate to "progressive** stall. 

Because i/^c(^) 0, Eq. (2.13) shows that A < 0, Thus, the trace 

of g in Fig. 2.2 is reversed in sign, and now it is the cell **leading 
edge** that is the steeper. The proper value of ri can be found from 
Eq. (2.19), if gjjj is fixed. 

Progressive Stall 

Pursuing this idea further, we can relate gnnn ^^d A$ from Eqs. 
(2.27), (2.28) to arrive at the criteria for limited, or progressive- 
type stall: 


S > ^ , a < 0 (2.29) 

i/l - (am/ n) ^ 

S > m , a > 0. (2.30) 
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These equations indicate that progressive stall is favored if 

1. S is large. Experience confirms that a delayed ’’break,'* 
and hence large S favors progressive stall. 

2. m is small. That is, an abrupt expansion at exit as opposed 
to a long, straight discharge would favor progressive stall. It is not 
known whether experience confirms this prediction. 

3. n is large, if a < 0 (clockwise hysteresis). Large n denotes 
many stall cells, and experience is quite definite that progressive stall 
is associated with multiple stall cells (^) . We find that n is relevant 
only if a < 0 and, therefore, f > 1/2 (Eq. (2.22)). It is also usually 
found that progressive stall has this feature as well (see Ref. 2 or 
Table 1.1), and we conclude that clockwise hysteresis (a < 0, T < 0) 

is quite typical in these "soft-stall" cases. Theoretically, it has 
been argued ( 11 ) that lift hysteresis of a single airfoil is of this 
sense if its stall process is a progressive advance of trailing-edge 
separation as opposed to the more abrupt leading-edge stall. In the 
former case, lift curves are rounded, while in the latter case, an 
abrupt break occurs. Thus, there is some basis for expecting T < 0 
for progressive compressor stall. 

4. a is small, if a < 0. Actually, am and n should be considered 
together; Eq. (2.17) provides 


am _ Dm 
n “ 2UNnT ’ 


(2.31) 


This grouping suggests that large T favors progressive stall, as does 
a large product Nn. Perhaps it may be inferred that large n is especial- 
ly required for single-stage (N = 1) progressive stall, whereas multi- 
stage machines can have progressive stall which does not split up into 
multiple cells because Nn is large even if n = 1. Experimental results 
tend to support this idea. 

Fig. 2.6 shows velocity and pressure-rise fluctuations calculated 
for a limited oscillation of the sort just discussed; with a negative 
and ipc quite large and negative, n was taken to be 4, and the stall 
point M was considered to occur at $ = 0.8. The amplitude is only 0.053, 
but n is quite large because is large, and the traces are quite far 
from harmonic shape. Propagation speed is greater than 1/2, and axial 
velocity changes most rapidly at the cell leading edge, in contrast 
to the case when > 0. This result is in contrast to the known ten- 
dency of deep-stall wave forms to be steeper at the cell trailing edge 
(^) . Probably, the differing curvatures of ipc are important in this 
regard. Perhaps it would be useful to investigate this aspect of pro- 
gressive stall experimentally. 
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Limitations of the Analysis 


Multi-cell progressive stall is usually also part-span stall (_2) . 
The present analysis assumes two-dimensional flow; in effect, full- 
span disturbances are assumed. The importance of this inconsistency 
with observations is difficult to assess. 

The foregoing analysis does not show that rotating-stall oscilla- 
tions must extend to the allowed amplitude limits, or that the oscilla- 
tion will choose, so to speak, to adopt a zone number appropriate to 
progressive stall. Neither is it shown that the oscillation will 
become of the large— amplitude type just because it is not limited to 
remain above point B. These questions are ones for future research, 
probably to be answered by studying limit cycles emerging from Eq. 

(2.4). The present analysis is intended to suggest features to look 
for, in those more complete and inevitably more complex and confusing 
solutions . 

Finally, we should recall the criterion — 0 mentioned 

following Eq. (2.3). Of course, a throttle setting need not correspond 
to $'($) = 0. However, a purely period wave cannot occur, governed 
by a purely parabolic characteristic unless that is true, even if the 
motion is nonlinear. The reason for this result, which will be made 
clear in Part III, is essentially that mentioned in Part I; namely, 
that ip^($) would be an amplification coefficient in Eq. (2.3). If 

^ then any oscillation would grow until it encountered new 

features of the function beyond the range of the assumed para- 
bola, such as a resistance to reversed flow, and this encounter would 

lead finally to a limit cycle. If i|Jc(^) < 0, then the oscillation would 
presumably disappear by damping. 

Thus, if $ is intermediate between points B and M, rather than at 
M, an oscillation would presumably amplify to include some Involvement 
with the steep drop at B. This case is beyond the scope of this paper. 

CONCLUDING REMARKS 


Finite-amplitude rotating stall has been analyzed for cases in 
which it may be assumed that the constant-speed characteristic is 
parabolic in flow coefficient with rotating stall absent. It was found 
that a periodic solution (neither damped nor amplified) requires a 
horizontal slope of that characteristic at the given average flow. 
Presumably, if the slope were negative, the oscillation would be damped, 
but if it were positive, the oscillation would grow until new features 
of the characteristic curve (reverse flow or stall recovery) become 
influential. 

In "deep stall," if the characteristic is flat, the solution 
differs only slightly from the linear one, perhaps explaining why 
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propagation speeds from Part I agree so well with experiment. In 
this case (Figs. 2.4 and 2.5), the oscillation would presumably be 
limited by reverse-flow resistance, provided flow coefficient is quite 
low (less than half the coefficient at stall) . At larger flow coef- 
ficients, even if the characteristic is locally flat, interaction with 
the unstalled part of the characteristic will occur, and the resulting 
motion will become non— harmonic on that account. This is the subject 
of Part III of this report. 

The characteristic function ^^^.((j)) we have used is artificial in 
the sense that it is defined in the absence of rotating stall, and 
therefore is not normally measured in compressor tests. However, it is 
perfectly well defined physically, and could be measured in experiments 
designed to suppress rotating stall. It is defined for the compressor 
alone, and therefore one is free to manipulate exit and entrance 
conditions in order to prevent rotating stall, and also to eliminate any 
irrelevant hysteresis effects of a control throttle. 

In incipient stall conditions, the stall oscillation may be 
limited in amplitude by a strong reverse curvature of the characteris- 
tic, permitting a progressive type of stall. If the characteristic 
"turns over," so that the break of the speed line is delayed substan- 
tially below maximum pressure rise, then this limited solution is 
favored. If the force lag is negative (opposite in sense to the effect 
of Inertia) then two notable results appear; The propagation speed 
exceeds 1/2 U, and the limited-amplitude solution is favored for multiple 
cells. Because experiments usually show these same features, the 
inference seems warranted that lag (t) is indeed negative in the incipi- 
ent or progressive stall region. 

In the deep-stall region, nonlinearity associated with concave 
curvature of the characteristic provides a steepening of the stall-zone 
trailing edge, as experiments show, provided lag (t) is positive. 
Probably, then, lag should be taken positive in deep stall, and negative 
in incipient stall. Experiments are needed to establish these blade- 
passage lags. 

Theory and experiment are both needed to describe the entrance and 
exit flows, especially for the relation between axial and transverse 
velocity disturbances. The present theory shows that rotating stall 
is sensitive to these conditions (lumped, in the present theory, in the 
parameter m) , and practical control means could well evolve from such 
a study. 

A new type of experiment is needed to find the compressor 
characteristic in the absence of rotating stall, and free of irrele- 
vant throttle hysteresis. Experiments have suggested that exhaust to 
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low pressure through a high-resistance screen may suppress rotating 
stall even at fully-stalled flow rates, and also that unwanted throttle 
effects can be eliminated by lowering pressure behind the throttle 
(see Ref. Fig. 19). Therefore, it seems feasible in principle to 
measure 

The present theory shows that the characteristics with and without 
rotating stall are different (6 represents that difference) and thus 
that rotating stall affects average compressor performance. Ultimately 
one can envisage a theory in which surge and rotating stall are 
integrated, rather than being considered separate and independent 
phenomena. One would then be able to deal with engine dynamics on 
a fundamental basis. 

A nonlinear, unsteady model of rotating stall has in fact been 
presented by Takata and Nagano ( 12) for single rows and pairs of rows; 
in it they permit the upstream potential flow to be unsteady, and make 
the required numerical calculations. Their assumptions about blade-row 
losses and pressure rise are quite different from those of the present 
analysis, and seem quite unsuitable for multistage compressors. Their 
chief concern was to analyse the nonlinear growth of disturbances into 
rotating stall of single and double rows. 

In the next and concluding part of this report, oscillations will 
be allowed to interact with the entire characteristic diagram, so that 
resistance to reverse flow and ’’unstalling” affect the disturbance 
pattern. 



Fig. 2.1. Constant-speed compressor characteristic in the absence of 
rotating stall. 
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Fig. 2.2. Nonlinear velocity disturbances, for various values of 
amplitude parameter n, from Eqs. (2.6a,b). 
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Fig. 2.3. Sketch of the increase of stream- tube area from A 2 to A 2 
upon entry into an IGV passage, owing to transverse dis- 
turbance velocity u. 




Fig. 2.4. Velocity and pressure disturbances for a flat characteristic 
= 0) in deep stall. Stall zone is 180 < 0 < 360 degr, 
where g < 0. 
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PART III - LIMIT CYCLES 


SUMMARY 


A theory of rotating stall, based on single parameters for blade- 
passage lag and external-flow lag and a given compressor characteris- 
tic yields limit cycles in velocity space. These limit cycles are 
governed by Lienard's equation with the characteristic playing the 
role of nonlinear damping function. Cyclic integrals of the solution 
determine stall propagation speed and the effect of rotating stall on 
average performance. Solutions with various line-segment characteris- 
tics and various throttle settings are found and discussed. 

There is generally a limiting flow coefficient beyond which no 
solution is possible; this probably represents stall recovery. This 
recovery point is independent of internal compressor lag, but does 
depend on external lags and on the height-to-width ratio of the diagram. 
Tall diagrams and small external lags (inlet and diffusor) favor 
recovery . 

Suggestions for future theoretical and experimental research are 
discussed. 


INTRODUCTION 


In Part I, it was shown that a simple time lag of blade-passage 
force, accumulated in a multistage compressor and suitably matched 
with lags arising in entrance and exit flows, permits rotating stall 
as a kind of eigensolution of the equations for weak disturbances. A 
prediction of stall propagation and speed results, and it agrees well 
with the experiment discussed by Day, Greitzer and Cumpsty (D.G.C.) (^) . 

In Part II, finite amplitude disturbances were considered in 
conjunction with compressor characteristics (upstream total to down- 
stream static pressure rise, vs axial flow coefficient) which have 
parabolic shape over the range of fluctuation. It was found that, if 
the characteristic diagram is flat, the linear, harmonic solution of 
Part I applies, presumably limited by resistance to reverse flow. If 
the diagram is concave the disturbance wave is distorted, with the 
reduced-flow (stall) zone being steeper at its trailing edge, as 
observed by D.G.C. If the diagram is convex, as at stall inception, 
the disturbance is of limited amplitude, especially for multiple-zone 
stalls when the lag is negative. In this case, propagation speed would 
exceed 50 percent of wheel speed. These last results may relate to 
"progressive stall" often observed near the inception point. 
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None of the foregoing results indicate "recovery"; clearly, a 
more complete diagram must be used, incorporating unstalled and probably, 
reversed-flow segments. Use of a more complete diagram should also 
provide estimates of the effect of the presence of rotating stall on 
the average compressor performance. This, then, is the theme of this 
part, to analyze the recovery phenomenon and other results which would 
come from using a complete performance diagram, and allowing oscillations 
to enter both unstalled and reverse-flow regions. We will find "limit 
cycles" which are more complex than the circles in phase space which 
arose in Parts I and II, and, under certain conditions, these limit 
cycles will appear to describe "relaxation oscillations", which we 
will argue represent recovery from rotating stall. Throughout, simple 
assumptions about internal and external lags will be maintained as in 
Parts I and II. 


THE CYCLE EQUATION 


We return to the basic equation (2.1) of Part II. It permits 
finite amplitude disturbances of axial velocity (g) and transverse 
velocity (h) at entrance, as functions of angle (0) around the com- 
pressor. Assuming zero entrance loss at the IGV entrance, 


Xg' (6) - mfh + 'F - = 0. 


(3.1) 


Lags are embodied in the parameters X and mf, where f is propagation 
speed coefficient, and m reflects inlet and outlet flow lags; if the 
exit flow is a sudden expansion, and entrance flow is potential, then 
n = 1. If the exit is a straight duct, then we know only that m = 2 for 
small disturbances as well, with m = 2 being appropriate for compressors 
followed by long diffusers. Research is needed to give a more complete 
picture of exit lag effects. Here, we will show the importance of m, 
or the phenomena it imperfectly represents, for rotating stall. The 
compressor lag is embodied in the parameter X, also related to f 
according to Eq. (3.16) to follow. 

The last two terms of Eq. (3.1) refer to the compressor diagram, 
sketched in Fig. 3.1. To define imagine a test in which c}> is 

varied from design to reverse flow, but rotating stall is suppressed 
by forcing axial velocity to remain circumferentially uniform (a high- 
drag screen at exit would work, in principle) . In that case, would 
be the characteristic measured in the absence of rotating stall. An 
operating point of the compressor (with axial velocity no longer 
constrained to be constant) would lie on a throttle line, but at 4^, 
either above or below the characteristic. Thus, 6 represents the effect 
of rotating stall on performance at a given average flow coefficient 


41 



4> = <E». It will be convenient to define the characteristic function 
relative to the intersection of ({)-$. Thus, 


¥ - iP^(O) E 6; E G(g) 

and Eq, (.3.1) becomes 


(3.2) 


Ag’(O) ■■ mfh - G + 6 = 0. 


(3.3) 


This is the equation to be solved in this part. 

If one finds the disturbance functions g(0) and h(0), the dis- 
turbance of static-pressure rise through the compressor can be found 
from Eq. (2.18) . Aside from constants, relative to pU^ it is 

Cp = -mfh + ^>g + (g^+h^) . (3.4) 


This formula includes pressure recovery at the IGV entrance. 
Relation of h to g 


In order to solve Eq. (3.3) it is necessary to connect the func- 
tions h and g. If the upstream flow is potential flow in a straight 
duct, then h(0) is the Hilbert Transform of g(0), given for periodic 
functions by Eq. (1.17): 

I A— 

h--^/ g'(5)ln|sln ^ I d?. (3.5) 

II 0 ^ 

Because g is defined as a difference from average axial velocity, the 
average i5gd0 =0. It can be shown from Eq. (3,5) that it is also true 
that the average tfhd0 = 0. Physically, this means that no swirl can 
be created in the entering potential flow. 

Although Eq. (3.5) could be numerically incorporated into Eq. (3.3) 
(as Takata and Nagano ( 12 ) did for their model) , it would be desirable to 
use some simpler relation, such as the one which is valid for single-cell 
stall and small perturbations, according to Eq. (1.19), with n = 1: 

h’(0)=-g. (3.6) 
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This equation is true for the first terms of Fourier series derivable 
from Eq. (3.5); in other words, it is exact if g is purely harmonic. 
Because we are concerned with "deep stall" in this report, for which 
experiments (3,14) show single-cell stall, we will assume n = 1. 
Actually, results for higher modes can be inferred by substituting the 
appropriate multiple of wheel diameter in the parameter a of Eq. (3.16) 
to follow. 

The nonlinear solution, Eq. (2.6b), can be manipulated to yield 


h' - - (1 - 2 a h)g 


(3.7) 


and one could consider making more general use of such a nonlinear 
relation in Eq. (3.3), rather than Eq. (3.6). 

We shall see that a "blockage" function has particular relevance 
to recovery. That is, g(6) might be a step function as shown in Fig. 

3.2. The example chosen has a "blockage" of 1-3 20 percent. Such 

a "relaxation oscillation" will provide a severe test of the reasonable- 
ness of Eq . (3.6). Substitution into Eq. (3.5) gives the corresponding 
h(e): 

sin(y - tt3) 

with A being the amplitude of the g-oscillation. Eq. (3.8) is shown 
on Fig. 3.2; the result is singular at the two values of 0 where g 
jumps. For comparison, we show Eq. (3.6). Apparently, Eq. (3.6) 
represents the correct Hilbert Transform quite well on the average. 

Eq. (3.7) would be inappropriate because it would not have the correct 
anti- symmetry about the midpoints of the stalled and unstalled zones. 

We are therefore emboldened to adopt Eq. (3.6) for general use in Eq. 
(3.3); it is correct in the small perturbation (harmonic) limit, and 
has the correct anti-s 3 mnmetry as well as a reasonable amplitude in the 
"blockage" limit. 

A more complete justification of the use of Eq. (3.6), and a simple 
improvement which can be applied, is given in Appendix A. 

L ienard’s Equation 

Eq. (3.6) permits the transformation of Eq. (3.3) into a generali- 
zation of the classical van der Pol equation ( 10) , with the character- 
istic function G(g) playing the role of damping function. From Eq. (3.6), 


(3.8) 


h = 


7t3(1-3) 


In 
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g'O) = - g dg/dh. 


(3.9) 


In phase space, Eq. (3.3) becomes 


Xgg ' (h) + mf h + G - 6 = 0 , 


The following transformations and shift or origin (and sign) of h 


_ G(g) ■ 


; 6 E 


u = Jlm ^ = A- /T 

mf » S mf / mf ^ 


yield a new expression for angle 0, 


d0 = - 


and yield the Lienard Equation ( 10) 


^ V /\ /% , /S. 

gg (h) + h = G(g) . 


Eq. (3.12) provides, for a complete cycle. 


A 

^ = 2tt / ^ 


and the requirement c^hdG = 0, where Eqs. (3.11) and (3.12) are used, 
becomes 

<P = (27T / "Y") 6 . (3 


If G were given, Eq. (3.13)- could be solved, and mf/X found from 
Eq. (3.14). However, A and f are related (Eq. (2.7), neglecting effect 
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of guide vanes for present purposes) by the definition 


X = - - f); a E D/(2 UNt) (3.16) 

a L 

where a is a reciprocal compressor time lag, in effect (Eq. (2.16)). 
Therefore, f itself can be found. Then, the characteristic displacement 
6 (or 6) comes from Eq. (3.15). Actually, of course, G(g) would normally 
be given, rather than G(g), and therefore Eqs. (3.13) and (3.14) 
would have to be solved simultaneously. 

Under the foregoing transformations, Eq. (3.4) becomes 

Cp = (mf)^(^i-S) + /x ^ ^ (h-^)^]. (3.17) 


Fig, 3.3 sketches the usual geometrical construction for solving 
Lienard's equation (10) . The characteristic is reoriented (compare 
Fig. 3.1) to correspond to the classical damping function. Hereafter, 
in discussing solutions of Eq. (3.13), we will use this convention 
without further comment. At a point © during an oscillation, where 
the values and g apply, one locates point , so that the line (^ - 
(§) has a length fi-S(g) . The line d) - (^ has the length g. The 
negative reciprocal of the ratio of these lengths is, by Eq. (3.13), 
just g' (li) or the slope of the trajectory of the oscillation. Thus, 
that trajectory is locally perpendicular to the line (^ - (^ , as 
sketched in Fig. 3.3. This geometrical method gives a quick sense of 
the progress of any oscillation in phase space. 

Line-Segment Characteristic 

In the rest of this paper, we will consider the function G to be 
various special cases of the line-segment diagram shown in Fig. 3.4a. 

The quantities IB - IF are parameters of the diagram to be specified. 

The throttle setting is specified by 3. We would usually consider the 
origin of <j) to be at 3 = 0, so 3B = 

A particularly simple but nontrivial case is sketched in Fig. 3.4b. 
It is a ’’vertical" diagram in which the unstalled leg is vertical 
(D = 0) as is the reverse-flow leg (E = 0) . The deep-stall leg is 
horizontal (C = 0) . The circular arcs shown dashed in the sketch will be 
discussed later. 
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SOME LIMIT-CYCLE EXAMPLES 


For purposes of illustration, we first pose a series of ''indirect" 
problems, in which we specify G(g), find the resulting oscillation, 
determine f, and then infer the function G(g). We note that Eqs. 
(3.13-15) are invariant to a uniform change of scale. Thus, for example, 
we may imagine our diagram of G to be scaled so that its width is 
exactly 1; in other words, so that g = g/B. On this basis, we assume 
our diagrams of G to have width 1, with the scale of h suitably adjusted. 
Fig. 3.5 shows a particular choice of diagram shape, double carets 
signifying the scale adjustment just discussed. 

1. Suppose 3 = 0.5, and an oscillation begins at a point such as 
P. An analytic solution of Eq. (3.13) can easily be found for any 
part of an oscillation for which G(g) is a straight line, and these 
solutions can be pieced together to form a continuous trajectory in 
phase space. In this case, the trajectory is at first a growing spiral 
(the construction of Fig. 3.3 suggests this result directly) which grows 
until it encounters reverse-flow resistance, and then becomes momentarily 
"unstalled." After about three circuits, the oscillation becomes self- 
duplicating; that is, a "limit cycle" of the sort familiar in nonlinear 
mechanics, and in fact which is classically Illustrated by Lienard’s 
equation with an S-shaped damping function. The usual compressor 
diagram, when properly oriented, is of that classical shape. 

The initial amplification in this example is due to the deep- 
stall slope being positive. If it were negative, the oscillation would 
locally damp. In any case, the peak of G being to the right of the 
reverse-flow corner is amplifying, for an oscillation which is already 
strong enough. The slopes of unstalled and reverse-flow legs being 
negative and (locally) limiting. Thus, the oscillation, as it grows, 
ultimately settles into a limit cycle. 

2. The limit cycle depends on the throttle setting ( 3 ) for a 
given diagram. Fig. 3.6 shows the effect. Various symbols on the cycle 
traces serve to relate each trace to the value of 3 indicated in the 
legend and along the deep-stall line where each cycle is centered. It 
is noteworthy that if 3 is somewhat less than 0.5, the cycle is "caught" 
by the reverse-flow leg- As 3 approaches 1, the "release point" moves 
closer to the maximum-G corner, and, in fact, the unstalled part of the 
cycle tends to hug the unstalled leg of the diagram. If 3 is slightly 
greater than 1, the unstalled leg actually catches the trajectory before 
it reaches the corner, and a rapidly damped oscillation follows (because 
the unstalled leg has negative slope) . 

Table 3.1 shows the results when Eqs. (3.14) and (3.15) are used 
to find f and hence G(g) and 6. We assume m = 2, and a = 1, in order to 
make the example definite. The actual F is only about 43 percent of B; 
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rather a lightly-loaded compressor. The value of f is about 0.26, 
tending to increase as 3 increases. As is obvious from Fig. 3.6, the 6 
needed to bring the average h to zero is greatest when 3 = 0.5; that 
cycle lies more to the right than the others. 

Fig. 3.7 shows the profiles of g(0) which may be derived from 
Fig. 3.6 and Eq. (3.12). The wave is propagating to the right. Thus, 
the fraction of circumference with reduced flow becomes less as 3 
increases, as one would expect. There is a tendency, not very strong, 
for the "trailing edge" of the reduced-flow (stall) zone to be steeper 
than the leading edge, for the lower flow rates, as observed by D.G.C. (^) . 

3. The limit cycle depends on the diagram "tallness." Fig. 3.8 
shows how the limit cycle changes when 3 is kept at 0.9, but the scale 
of G is changed by factors of 2 and 3 beyond the case shown in Fig. ^3. 6 
(and repeated in Fig. 3.8 for comparison). Obviously, as scale of G 
increases, the cycle follows the unstalled leg very closely indeed, 
releasing just at the top point of the diagram in the limit. The 
reverse-flow leg is also followed closely, and, again, release occurs 
close to the lower corner of the diagram. These effects of diagram 
tallness will be of great significance in our subsequent discussion 
of recovery. 

The last three rows of Table 3.1 show the completed calculations. 

Even in the most extreme case, IF is only about equal to IB. In that 
case, f has increased to 0.43, and 6 has increased dramatically, as 
Fig, 3.8 suggests it should. 

Fig. 3.9 shows g(0) for the most extreme case, to be compared 
with the profile 3 = 0.9 in Fig. 3.7. The stalled angle is reduced, 
and the expected steepness of the stall trailing edge is very definite. 

The intercepts of a horizontal line (G = 0) with the unstalled and 
reversed legs are shown (unstall g" and "reverse g") to indicate how 
well this profile conforms to blockage ideas, namely that the flow is 
either stalled or unstalled. The concept seems not to be very useful 
in this case. Also shown is the pressure-rise disturbance Cp, obtained 
from Eq. (3.17), with B taken to be 1. 

We proceed now to analysis of the direct problem, in which a 
simple "vertical" shape of G(g) is given, and we will pay particular 
attention to indications of stall recovery as the throttle is opened. 

IDEAL-DIAGRAM SOLUTIONS 


Vertical Diagram 

The "vertical" diagram of Fig. 3.4(b) has only F and B as parameters, 
and, because the equations are invariant to overall scale change, only 
the single parameter F/B is significant, along with throttle setting 3. 
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The limit cycle in this case is simply described in the fi,g plane. 
Reference may be made to Fig. 3.4(b), imagining it to be transformed 
to h and g, or as if mf = 1. Eq. (3.13) tells us that the cycle 
follows the unstalled leg (g’ =0) until the peak of F is reached 
(h = ^ay) . The cycle leaves that corner, and is governed by G(g) = 0. 
In that part of the cycle, Eq. (3.13) is satisfied by a circular arc 
centered at (0,0). It is shown dashed on Fig. 3.4(b). The arc will 
intersect the reverse-flow leg as shown, provided the diagram is 
tall enough. 


In the h,g plane, ’’tallness" will depend on f as well as F/B, 
because of the transformations, Eqs. (3.11) and (3.12). The height- 
to-width ratio will be 


F/B 



in the fi,g plane. Thus, if f is near 0.5, so that A is small, then 
any diagram will be "tall" in the h,g plane, and, incidentally, the 
cycle will always reach the reverse-flow line. 

The cycle will then follow the reverse flow line until the corner 
is reached at h = 0. Then another circular arc governs the return to 
the unstalled condition. This second arc will reach the unstalled 
line only if 3 >0.5. If 3 <0.5, the arc will simply complete a 
circle back to the starting point. This solution corresponds to the 
purely harmonic wave (limited by reverse-flow resistance) discussed 
in Part II. 

It is a simple matter to evaluate the cyclic integrals, Eqs. (3.14) 
and (3.15) in these cases. Along the arc segments, the differential 
dh/g is simply the angle increment along the arc. T heref ore, according 
to Eq. (3.12), the real angle change, d0 , is simply /A/raf times the 
angle of arc in the h,g plane. An immediate conclusion may be drawn: 
Fig. 3.4(b) shows that the arc from unstalled to reverse condition 
generally subtends a smaller angle than the return. Therefore, the drop 
into stall is more abrupt than the recovery. Since the stall zone is 
traveling in the positive 6 direction (Fig. 1.4), the drop into stall 
is the "trailing edge" of the stall zone, in D.G.C.’s terminology (_3) , 
and, in keeping with their observations, the trailing edge of the zone 
is predicted always to be steeper than the "leading edge." 

The foregoing simple calculations were made for various F/B, A, 
and mf, subject to the requirement of Eqs. (3.14), and finding 6 from 
Eq. (3.15). Typical results are shown for F/B = 1 on Fig. 3.10. The 
solid-line results do not depend on compressor lag (a); only when the 
definition of A (Eq. (3.16)) is introduced, overlaid on Fig. 3.10 as 
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dot-dashed lines, does a become a factor. The intersections of dashed 
and solid lines provide the solution shown in Fig. 3.11 for the special 
cases of m = 2 and a = 0.3 and 0.6. 

These solutions describe two regimes: The harmonic wave, for 

which 6=0, when 3 < 0.5, and a rather different wave when 3 > 0.5. 

The velocity, pressure-rise, and phase diagrams for these two cases 
(with a = 0.6) are shown as the first and second rows of Fig. 3.12. 

The steeper ’’trailing edge” is evident in g(0) for 3 > 0.5. 

It is obvious from Fig. 3.10 that the external lag parameter (m) 
is more significant than the internal one, especially for the larger 
values of g. The speed f is restricted to be less than 0.5. Therefore, 
if ra is small, say 1, rotating-stall solutions might be forbidden 
entirely. Even if m = 2, no solution is possible when 3 > 0.801. In 
the next paragraphs, this point will be discussed more fully. 

The Recovery Limit 

When m = 2, a limit is reached when 3 = 0.801. At that limit, 
which is independent of a, the speed f goes to 0.5, and so does 6/F. 

The third row of Fig. 3.12 shows the corresponding traces. The axial 
velocity becomes a pure step wave and the circumferential velocity 
is a saw tooth. (This type of "relaxation oscillation” was discussed 
earlier in connection with Fig. 3.2.) 

The phase diagram becomes a rectangle because X has approached 
zero. From Eq. (3.10) one would expect dg/dh to become infinite when 
X ->■ 0. Another way to see this result is to recall that the cycle 
includes arcs in the h,g plane, in which the diagram is infinitely 
"tall” when X ->■ 0. That is, while g changes by B, Eq. (3.12) says 
that g changes only to order B/X when X is small. Thus, the arcs depart 
to a vanishing degree from = IF/(l/2)m when stalling or h = 0 when 
unstalling. The result in the original frame of reference is a rec- 
tangular cycle. 

Pursuing further the effort to interpret the limiting solution 
physically, we may focus attention on Eq. (3.1A). Essentially, it is 
the requirement that the cycle complete itself in the angular distance 
2 tt. Now if we specify a large 3» then g is a small constant quantity 
on the unstalled leg of the cycle. The coordinate fi must change by a 
certain amount essentially fixed by the height of the diagram (F) . Thus, 
the contribution of the unstalled leg to the integral g5dfi/g will be 
of the order F/ [ (1-3)B/X] , and Eq. (3.14) suggests that if 3 is too 
near 1, the unstalled leg will use up all the allowed range of 0. 

Clearly, 3 he permitted to be largest when the other terms of the 
integral are minimized. This occurs when the stalling and unstalling 
arcs have the^least angle in the fi,g plane, namely zero when X -> 0, and 
the negative-h portion of the unstalled leg is minimized; again, when 
X 0. In the X -► 0 limit, the reverse flow leg is longest, but g is 
larger there, and so the contribution to 3 is less important. 
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Appendix B provides further physical discussion of the relaxation- 
cycle limit as a matter of "detailed balancing" of pressure changes in 
a stage. 

In summary, as the throttle is opened, a limiting condition is 
reached such that the cycle spends too much time unstalled to permit 
the required pressure changes of stalling and unstalling to occur. 

This description would seem to fit the phenomenon of "stall recovery" 

O ) . It is significant that 6 is quite large (about (1/2)F) when the 
critical 3 ~ 3cr is reached. If this is true, a discontinuity of V is 
predicted to occur when 3 = 3cr* geometry of Fig. 3.1 suggests that 

a slightly more open throttle line than that for 3cr would have no nearby 
performance line to intersect, and a jump across to the unstalled line 
(i.e., recovery from stall) would have to occur. 

The value 3cr ~ 0,801 applies only for a vertical characteristic 
with F/B = 1 and m = 2. The general results for a vertical characteristic 
are 




iF/B 

Trm 



(3.18) 


which appear as the lowest curve on Fig. 3.14. If the diagram is less 
"tall," 3j,j. can approach 1. On the other hand, if F/(B exceeds (l/4)7rm, 
no rotating stall other than the harmonic variety can occur. 

In general, a tall diagram (large F/B) and a small value of m 
favor recovery. The role of m is particularly significant. We have 
inferred that m is about 2 for a compressor with a straight outlet. 

If m were reduced to 1 by using a sudden expansion, rotating stall 
could be largely suppressed, according to the present analysis, when 
F/B is greater than about 0.8. Thus, changes in external lag mechanisms 
could be very effective for relatively "tall" diagrams, it would seem. 

The Step Diagram 


Another sort of ideal, simple diagram would continue to have C = E = 
0, but D = <» rather than 0. Fig. 3.13(a) shows such a "step" diagram, 
together with a phase diagram in the h,g plane (as if mf = 1). Such 
an extended unstalled plateau seems a possible model for lura and 
Rannie’s diagram (13) , but is not typical of modern compressors. 

The results show a lower, rather than an upper limit for rotating 
stall. It is 


3 


cr 


2 F/B 

TTm 


(3.19) 
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shown in Fig. 3.13(b). If the diagram is tall enough, there is a 
forbidden gap between 3 = 0.5 and 3d- where a rotating-stall solution 
is impossible. This is true only if F/B exceeds (1/4) inn, which is 
quite tall. For smaller values of F/B, there is no critical case, and 
solutions for f are quite far from 0.5. For example, if F/B = 0.5, m = 2, 
and a = 0.6, the speed f would remain in quite a narrow range; 0.35 
at 3 = 0.5 and 0.25 at 3 = 1.5. 

RECOVERY FOR LINE-SEGMENT CHARACTERISTIC 


We now return to the more general line-segment characteristic 
shown in Fig. 3.4(a), and look for the effects of nonzero C, D, and E. 
The general solution is easy to find, but hard to discuss because there 
are so many parameters to consider; the previous section describing 
Figs. 3.6-9 will suffice for illustration. If we consider only 
recovery, simplifications are available; f approaches 0.5, X approaches 
0, and compressor lag (a) becomes irrelevant. 

As X 0, we have seen how the transformation Eq. (3.12) makes 
the diagram "taller” in the h,g plane without limit, so that, as in the 
vertical diagram case, the installing and stalling legs are lines of 
constant fi. The deep-stall slope C becomes irrelevant except in 
finding 6, provided C is positive. If C is negative, the cycle would 
presumably follow the deep-stall leg and release at the corner g = 
(1-3)(B. We do not consider this possibility any further in the present 
paper. With (C > 0, the only change from the vertical-diagram analysis 
is that the integral of dfi/g must be carried out along the slanted 
lines defined by D and [E. When this is done, the following formulas 
for recovery result: 


TTm 1-3 3 

cr cr 


(3.20) 


^_cr ^ ^cr 

^ (D/B) ^ 


1 / -I I D/B X IB 
cr 


+ ind + - |l . 

(E/B)^ ^cr ^ 


(3.21) 


These formulas reduce to Eq. (3.18) when E and D are set equal to zero. 

The formula for 6 is valid even if (C is not zero, so long as 6 is measured 
from the reverse-flow corner at g = -3B. 
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Fig. 3.14 shows the effect of "slant," choosing D = E for con- 
venience. Clearly, slant increases the value of 3crJ other words, it 
delays recovery. It remains true, however, that a tall diagram with 
small m will favor recovery. Fig. 3.14 shows how the maximum value of 
F/mJB, which occurs at 3cr “ 0.5, depends on the slant parameter D/B. 

The more slanted the characteristic, the taller may be the diagram 
before the rotating-stall limit descends to 3cr " 0.5. Fig. 3.14 also 
suggests that for realistic slants of about 0.5, recovery occurs in a 
fairly narrow range of 3 (between 0.9 and 1.0 for F/B < 0.7m). In 
D.G.C. (3), it is suggested that experiments indicate recovery when a 
certain critical blockage (about 0.7) occurs. The present analysis is 
not quite consistent with this idea because 3cr same as 

blocl^ge. If 3cr near 1, and D/B = 0.4, then the two ideas do 
happen to agree. 

The dashed lines of Fig. 3.14 show that 6^^ is not greatly 
affected by slant unless 3cr quite near 1. Of course, that is the 
region where recovery is generally predicted. 

Further results of Eqs. (3.20,21) are provided in Appendix C, 
together with a set of recovery predictions which emphasize that only 
the top and bottom points of the diagram are needed to fix the relaxation- 
cycle limit. 


COMPARISON WITH EXPERIMENT 


Comparison of the present theory with experiment is handicapped 
by ignorance of the rotating-stall-free characteristic i|J(,(^)* In 
principle, it could be measured, but heretofore there has not been 
sufficient reason to do so. Nevertheless, in Part I it was shown that 
theoretical predictions of propagation speed (f) agreed well with 
experiment. Here, we are concerned chiefly with recovery, and we may 
appeal to D.G.C. (^) for evidence. 

Recovery 


Fig. 3.15 shows characteristics measured and reported by D.G.C. 
These appear as solid lines, with the recovery point Indicated by a 
circle. Three different loadings are included, all at 50 percent 
reaction. For theoretical predictions, we have assumed that the shut- 
off head describes a constant value of deep-stall iJJq, and calculated 
values of 6 are applied relative to that level of 

Therefore, if vertical diagrams were used, the theory would follow 
the dashed lines shown in each of the three cases. The average slant of 
the unstalled leg gives a basis for estimating the recovery point by 
use of Eqs. (3.20) and (3.21). These points are triangles, and are pre- 
sumably better estimates than the diamonds which terminate the vertical- 
characteristic results. 
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The "low-(J)*" (Case I) results are not very Informative. The 
"intemediate-cj)*’' (Case II) results show excellent agreement, even 
including a suggestion of transition from harmonic to limit-cycle waves 
The ”high-(j)*" (Case IV) is reasonable as regards 6cr> quite poor 

as regards vertical shift of ' Probably, the slope C of the deep- 
stall should be taken as negative. If so, the integral j^dh/g should 
be recalculated. It seems doubtful, however, that such a revision 
would give such a low value of 'F at recovery. Perhaps the throttle 
plays a special role in recovery when the characteristic descends to 
such a low limit. 

In all cases, it is worth noting that recovery, both in theory 
and experiment, is close to the breakpoint of the characteristic. That 
is, 3 ~ 0.95 at recovery. Perhaps this may prove a more basic concept 
than blockage. 

Propagation Speed 

Fig. 3.16 shows the propagation speeds which correspond to the 
Cases I, II, and IV of Fig. 3.15. These speeds, given in Day's Thesis 
(14) , relate to theory just as recovery does in Fig. 3.15, if we recall 
that f should approach 0.5 (A 0) at recovery, according to Fig. 3.11. 

The low-(f)* Case I does not show the predicted rise of f toward 0.5 
but the measurements do not extend to recovery; perhaps the rise would 
indeed occur if recovery were more closely approached. 

The intermediate-(})* Case II does show the predicted rise; perhaps 
the good agreement in this case is related to the good agreement of 
recovery prediction and experiment shown for Case II in Fig. 3.15. The 
high-<j)* Case IV shows the least agreement; the experiments show no rise 
all the way up to recovery. 

lura and Rannie ( 13 ) also report a propagation speed which remains 
nearly constant right up to recovery. In their experiments, the com- 
pressor has an extended plateau of pressure beyond the stall point, 
somewhat like the steep characteristic of Fig. 13. In that case, we 
would in fact predict a nearly flat value of f. 

Further comparisons of theory with experiment will require, on the 
one hand, a more thorough theoretical exploration of the parameters 
of the problem, and on the other hand, new measurements of rotating- 
stall-free compressor characteristic maps. 
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CONCLUDING REMARKS 


The problem of rotating stall has been difficult largely because 
it occurs in such a complex parameter space. Perhaps the best service 
theory can provide in such a situation is to suggest priorities among 
these parameters, rather than to predict the resulting phenomenon in 
detail. The present set of three reports is set forth in this spirit. 
The most important results seem to be the following: 

1. The recovery from rotating stall seems basically predictable. 

It depends on the compressor's steady characteristic; primarily its 
"tallness" (the height of the drop into stall compared with the diagram 
width), and secondarily the slope of the unstalled characteristic. 
Tallness favors early recovery. The slope of the reverse-flow resis- 
tance is of minor significance, as is the slope of the deep stall line, 
provided that slope is positive. If it is negative, the throttle 
curve may be important; otherwise, the throttle characteristic is 
unimportant. 

Recovery does not depend on compressor lag (dynamic hysteresis) , 
but it depends very much on the lags originating in the entrance and 
exit flows. Small external lags favor recovery. It seems likely 
that the most directly effective design change would be to reduce the 
lag in the exit duct. Appendix D provides a brief analysis which 
indicates that a converging entrance would provide a reduced external 
lag (smaller value of m) also. 

2. The characteristic with rotating stall suppressed must be 
found experimentally to form a proper basis for further theoretical 
refinements. Discharge to a very low pressure through a dense screen 
at exit would force circumferentially uniform axial flow, and suppress 
rotating stall. Also, it would be important to suppress any effects 
of finite throttle-line slope, such as steady hysteresis about stall. 
Again, discharge to a low pressure would tend to steepen the throttle 
lines in the necessary way. 

3. Research is needed on entrance flow fields and, especially, 
diffusers with transversely non-uniform flow. Dynamic (lag) effects are 
very important, more so than compressor blade-row d 5 uiamics. 

4 . Theoretical limit-cycle analysis can easily be extended to 
include a wider class of hypothetical characteristic diagrams, and a 
certain "exact" entrance and exit-flow relation. It is felt that these 
should proceed in parallel with the studies discussed in 2 and 3 above - 

Appendix E contains more detailed discussions of the implications 
of the theory of this report for future experimental, theoretical and 
computational research. 
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TABLE 3.1. RESULTS FOR LIMIT-CYCLE EXAMPLES 


Shapes of G(g) » as in Figs. 3.5 and 3.8, Refined by reference to^ 
Fig. 3.4: fi = 1; 6 = = 9/16; ^ = 0.9; C = 0.2. In Fig. 3.8, F 
and £ multiplied by 1, 2, 3. m = 2, a = 1. 


Symbol 

3 

/ X 

§ 

f 

/Xmf 

F/B 

6/F 

A 

0.25 

1.028 

0.043 

0.26 

0.50 

0.45 

0.048 

0 

0.50 

1.076 

0.221 

0.26 

0.60 

0.45 

0.246 

□ 

0.75 

1.084 

0.190 

0.26 

0.50 

0.45 

0.121 

O 

0.90 

1.151 

-0.032 

0.29 

0.50 

0.45 

-0.036 

V 

1.00 

1.336 

0.068 

0.32 

0.48 

0.43 

0.076 

o 

0.90 

1.151 

-0.032 

Pi 

0.50 

0.45 

-0.036 

0 

0.90 

1.178 

0.387 


0.44 

0.78 

0.215 

o 

0.90 

2.400 

0.766 

0.426 

0.36 

0.96 

0.284 
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Fig. 3.1. Sketch of notation on compressor performance diagram, of 

total-to-static pressure rise flow coefficient (cj)) 
in absence of rotating stall. Operating point in rotating 
stall is 6 above steady diagram. 
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Fig. 3.2. Comparison of correct determination of h from g (Eq. (3.8)) 

with generally approximate Eq. (3.6), for a step wave ($ = 0.8) 





Fig. 3.3. Geometrical solution of the Lienard Equation (Eq. (3.13)). 



(a)General line-segment (b)"Vertical"diagram 
diagram (C=D = Z=0) 


Fig. 3.4. Line-segment characteristic diagrams, of function G(g) , 
oriented as in Fig. 3.3. 
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Fig. 3.5. Development of limit cycle, for 3 = 0.5. Oscillation begins 
at point P. After three cycles, negligible changes. See 
Table 3.1 for diagram parameters corresponding to assumed 6(g). 



Fig. 3.6. Effect of throttle setting (3) on limit cycles; for the 

same given 6(g). See Table 3.1 for the corresponding diagram 
parameters. 
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Fig. 3.9. Axial velocity (g) , static pressure-rise (Cp) disturbances, 
for 3 = 0*9 and the most extreme case of Fig. 3.8 (Symbol 0? 
see Table 3.1) . 



Fig. 3.10. X and 5 (shown dashed) vs mf for various B and IF/B = 1. 

The definition of X = (1/a) [ (l/2)-f ] overlaid (dot-dashed). 
Point P' gives the solution for 3 = 0.5 when a = 0.6, m = 2. 
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Fig. 3.11. Propagation speed (f) and performance rise due to rotating 
stall (6) shown dashed vs throttle setting (B) > for m = 2 
and F/B =1. No solution for m = 1. Note that B^r “ 0.801 
is independent of lag (a), and that 6=0 for B < 0.5. 
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Fig. 3.12. Theoretical velocity and pressure "traces'* and phase diagrams, for a particular vertical 
characteristic, and three values of 3 (Bqj- = 0.801), m = 2, a = 0.6. Arrow denotes 
direction of propagation. 
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Fig. 3 . 14 . Effect of slant characteristic on recovery throttle setting 
(3cr) and performance rise at recovery (6„r), shown dashed, 
for E = D. 
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II 
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Fig. 3.15. Comparison of theoretical predictions with cases I, II, and IV (50% reac 
Greitzer and Cumpsty (D.G.C.), Ref. 3. 
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APPENDIX A 


FURTHER ANALYSIS OF THE h-g RELATION 


In Part III, the relation h' = -g was used (Eq. (3.6)) in place of 
the exact potential flow equations for a straight entrance channel. That 
exact relation is given by Eq. (3.5) or the equivalent Fourier Series; 


If 


then 


00 



sin n6 + b 

n 


h = I 0 


cos n9 - b 

n 


cos n0) N 
sin n0 ) ^ . 


(Al) 


In Part III, Fig. 3.2 compares h' = -g and Eq. (3.5) for a step-function 
form of g(6). While the result is quite good generally, the simpler 
form fails near the jumps. Of course, if g(0) is harmonic, then h' = -g 
is exact. 


More recently, an improved relation was studied: 


h' = -Ag + Bg”. 


(A2) 


In effect, Part III assumed A = 1, B = 0. Now we ask whether there 
might be better choices of A and B. Substituting Eqs. (Al) , we find 


h' = - y (n) (a sin n0 + b cos n0) 
“ n n 


-Ag + Bg" = - y (A + Bn ) (a^ sin n0 + b^ cos n0) 


(A3) 


Apparently, the harmonic terms of these two series are identical, and 
we can match two sets of coefficients, satisfying Eq. (A2) to that 
extent . 


If only first and second harmonics are matched (n = 1, 2), then 
Eqs. (A3) give 


A = 2/3; B = 1/3 


(A4) 


while if only first and 


third are matched. 
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A = 3/4; B = 1/4. 


(A5) 


Thus, we infer that the best choices of A would typically be a bit 
less than 1, and B should be quite small. 


Recalling Fig. 3.2, we see that a choice like Eq. (A5) for the 
step-wave would be equivalent to lowering the slope of h(0) and adding 
a 6(6) function at the jumps; this would be a plausible improvement. 


Part II provides certain exact solutions of Eq. (Al) . 
(Eqs. (2.6a,b)) 


i0 


’’o " 


hi + igl E 


1 + ne 


i6 


. i0 
xe 


1 + ne 


i0 


They are 
(A6) 

(A7) 


where n is an arbitrary parameter denoting departure from harmonic shape. 
Fig. Al shows the exact values of hQ and hj^ gotten by differen- 
tiating Eqs. (A6,7), for a particular choice, n = 0.4. 

The exact solutions are contrasted with applications of Eq. (A2) for 
certain A, B combinations. It appears that A = 0.75, B = 0.25 is best. 
These comparisons are repeated for a more extreme choice of n = 0.7, 
and are shown on Fig. A2. In this case, A = 0.85, B = 0.15 seems best. 

We conclude that Eq. (A2) should represent Eqs. (Al) quite accurately, 
with A chosen in the range 0.75-0.85 and B in the range 0.25-0.15. 

Of course, changing the h-g relation will modify Eq. (3.10). The 
required change is trivial, however. Define a new variable z, 


g = -z ' ; h E Az - Bz". (A8) 

Then, Eq. (3.3) becomes 

(A-Bmf )gg’ (z) + Amfz + G(g) -6=0 (A9) 

instead of Eq. (3.10), and 


d0 = -dz/g; {^zdz/g = 0. 


(AlO) 
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Thus, limit cycles in g,z space may be found just as they were found 
in g,h space in Part III, with X replaced by X - Bmf and mf replaced 
by Amf. For example, recovery will occur not at f = 1/2, but at a 
slightly lower value of f : 


f = 


1/2 

1 + Bma 


(All) 
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Fig. Al. 


Evaluation of h-g relations by use of exact solutions 
Eqs. (2.6), for moderate nonlinearity, T) = 0.4. 
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APPENDIX B 


DETAILED BALANCING AT RECOVERY 


In considering the "recovery" results, one is struck by the predic- 
tion (borne out be experiment) of sudden jumps of flow coefficient at 
stall— zone boundaries, and one may wonder how this is possible. The 
following argument seems to explain the matter. 

Beginning with Eq. (3.3), and assuming h* = -g, we may write 

ah" + mfh + G(-h’) +6=0 (Bl) 


where X ~ (2 UNt/D) ( 1/2 - f) represents the "mass" of an oscillator, 
mf the "spring constant," and G(-h’) + 6 the "dashpot." 

We know that recovery is associated with X = 0, or, we may say, 
vanishing oscillator "mass." Certainly, a massless oscillator can 
undergo a jump-like displacement such as we wish to explain. In our 
problem, vanishing mass evidently means either vanishing T or f equal 
to 1/2. The first possibility is to be rejected as physically untrue; 
there is always flow inertia in blade passages. 

What happens when f 1/2? Referring to Fig. 1.3, and Eq. (1.3), 
we see that (for the assumed 50% reaction) when f = 1/2, the acceleration 
(p is equal and opposite in rotor and stator blade passages. Specifi- 
cally, a pressure drop is required in the stator to accelerate the 
fluid there, and an equal pressure rise is required to decelerate fluid 
in the rotor. These pressure changes will each be (theoretically) 
infinite for a step-change of g(0) , but they will cancel for a stage. 
Thus, a kind of "detailed balancing" applies to inertial effects in 
each compressor stage. This result shows an important respect in 
which single-row and multistage rotating stall must be profoundly 
different. Cumpsty and Greitzer (^) give another expression of the 
foregoing idea. 

It would be interesting to monitor the row-by-row pressure changes 
experimentally, to verify the balancing process just described. 
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APPENDIX C 


FURTHER ’’RECOVERY" RESULTS 


A more complete version of Fig. 3.14 may be made by application 
of Eqs. (3.20,21). Fig. (Cl) shows the resvilts when D and E are not 
necessarily equal. Clearly, unstalled slope is more important than 
the reverse-flow slope over most of the range of 3* The corresponding 
results for 6 (Eq. (21)) give the same conclusion; they are shown in 
Fig. (C2) . Fig. (C3) shows certain profiles of g(9) at recovery. 

The recovery analysis is independent of the shape of the charac- 
teristic between its extrema. This fact is emphasized in Fig. C4, 
which shows the recovery points relative to the extrema in a few cases. 
The displacement is to be measured from the minimum point, in the 
recovery limit. This result suggests that the rotating-stall-free 
characteristic is not very important, if we are chiefly concerned with 
recovery itself rather than approach to recovery. 
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Fig. Cl. Recovery requirements, as they depend on slant features of 
characteristic . 



3 





1 . 0 - 



Fig. C3. Recovery profiles of axial velocity vs 0, as they depend on 
slant features of characteristic. 
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Fig. C4. Recovery point locations (circles points), depending on 
locations of maximum and minimum points, and slants, of 
characteristic, for two values of external lag coefficient (m) . 
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APPENDIX D 


WEDGE-LIKE ENTRANCE FLOW 


If the entrance duct is straight, and a harmonic wave of axial 
velocity occurs. 


V = c sin 0 


(Dl) 


then, according to Eq. (1.19), 


u = c cos 6 


(D2) 


and we may ask if a different shape 
would produce a different, perhaps 
smaller, value of u in response to 
an imposed v. If so, the important 
parameter m would be reduced, and 
rotating stall would be inhibited 
from occurring. 

To illustrate this idea, we 
consider a convergent, straight- 
sided entrance, as sketched. The 
shape is two dimensional, and the 
v velocity varies in the third 
dimension x, or 6, where x = (1/2)D0. 



We imagine the flow to be axisjnnmetric and irrotational , so that 
a velocity potential satisfies 


(b + — d) + 5 =0. 

rr r r xx 


(D3) 


A disturbance at the compressor entrance, where ar = H, is assumed: 


<j) = cT(i;)sin 0; 



(D4) 


whence a Bessel’s equation is found: 





- T = 0 
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which is satisfied by the function T = ^q(C) (vanishing at !^ = «») , 
defined by the series 


Kq = -(0.5772 + In |)(1 + i i + ... (D5) 


Thus, at the face, 


-V = ^ K^(Osin 0 

Kq(C)cos e / . 

If C is small there, Kq ~ -(1/C) and Kq ~ -(0.5772 - In 2/^). 
rOt = H (see the sketch), C/2 = (l/a)(H/D), and 

_ .2c iv , Q 

V ~ (— — ) sxn 0 
u ~ (In -| - 0.5772)cos 0. 


(D6) 


Since 


Identifying the coefficient of v as c. 


u c (In 2/C - 0.5772) cos 0 


(D7) 


which is to be compared with Eq. (D2) . With 


1 = 1 H 

2 a D 


(D8) 


we see that a large angle cx suggests a small induced u, hence a small 
value of m, and a reduced tendency to rotating stall. Specifically, 
if H/D = 1/10 and a = 60®, Eq. (D7) yields 


u = 0.338 c cos 0 


(D9) 


which is only one- third of the straight-duct result. 

We conclude that convergence of the entrance region could be an 
effective means of control of rotating stall. 
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APPENDIX E 


IMPLICATIONS OF THEORY FOR FUTURE RESEARCH 


Characteristic Free of Rotating Stall 

In this report, it has been assumed that a characteristic speed 
line is known, which is apart from the existence of rotating 

stall, and is free of any hysteresis due to throttle-line slope. Of 
course, this is not normally measured, because, in experiments, rotating 
stall does occur, and Parts II and III show that there is a substantial 
average effect of rotating stall to be expected when the amplitude is 
large. Also, throttle lines are parabolic, in practice. We note 
that surge theories or theories which combine surge and rotating 
stall would assume knowledge of 

Two questions arise; (1) How important is it to know 
(2) how might it be measured? On the first question. Part III, as 
more fully interpreted in Appendix A, shows that for recovery predic- 
tion, only the unstalled and reverse flow legs of are needed 
(assuming that minimum pressure rise is at shut-off) , and these can 
be found well enough from experiments in which rotating stall is allowed 
Of course, reverse flow is a rather confusing subject in any case. 
Fortunately, as the cited discussion shows, the details of reverse 
flow are of slight importance for recovery. 

Thus, in the deep-stall range of 0 is important only for predict 
ing the complete deep-stall rotating-stall performance. One may suspect 
that it will prove unimportant for the ultimate theory — engine- 
recovery prediction. 

The foregoing argument empha- 
sizes the input needed from experi- 
ment to make a recovery prediction. 

An experiment intended to check the 
prediction of recovery may be 
another matter, however. In the 
latter case, throttle related hys- 
teresis could spoil the comparison. 

That is, the throttle line might 
osculate the speed line before 
recovery is reached, as indicated in 
the sketch. Perhaps this is what 
happens in the cj)* = 1.0 case of 
Fig. 3.15. Thus, it may be more 
important to achieve a steep (ideally, 
vertical) throttle line than to 
suppress rotating stall. 
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How might be measured, free of (a) rotating stall or (b) throttle 
effect? To accomplish (a) , it may be enough to prevent tangential 
flow velocities in the entrance or exit regions. In principle, honey- 
comb or vane guidance will work, provided vanes extend all the way to 
the supply and receiver reservoirs. Probably, some compromise is 
possible in that regard. A single stage should suffice, provided it is 
truly representative, because the overall lag in the compressor, which 
depends on number of stages, is unimportant when rotating stall is 
to be suppressed anyway. 

Another approach might be to add an unstalled compressor in series, 
closely coupled to the test compressor. However, even the combined 
characteristic would probably still have a sudden drop where the test 
compressor stalls, and it seems likely that the combination would still 
show come rotating stall. For purposes of rotating-stall suppression, 
the fore and aft honeycomb would seem the most feasible. 

To accomblish (b) (elimination of throttle effect) , two approaches 
again seem possible. The most direct method would be to drop the 
receiver pressure (or raise the supply pressure) . That would move the 
origin of the throttle parabola downward on the diagram, thus steepening 
the throttle line where it cuts the characteristic at a given flow 
rate. The range of effect obtainable by this method may be limited. 
Rather than lowering the throttle origin, it might be better to raise 
the characteristic by adding a compressor in series, before the throttle. 
It would not need to be closely-coupled. It should have a high pressure 
ratio, and, if possible, a rising characteristic as flow is reduced. 

A centrifugal compressor would be a promising possibility. A difficulty 
would be that the add-on compressor would have to be at least as 
powerful as the test machine. To accomplish (b) alone (allowing rotating 
stall) would mean that the full number of stages would be needed, to 
get the rotating stall right. 

In summary, it would seem possible to find ip^(6) , suppressing 
rotating stall and throttle-related hysteresis, by adding axial guide 
vanes extending into slow-flow regions ahead and behind the test com- 
pressor, and by adding a second compressor in series betwen the test 
compressor and throttle. The latter is probably more crucial if the 
only concern is "recovery." 

Entrance/Exit Dynamic Relations 

In this report, some rather far-reaching assumptions were made 
about inlet and exit flows, even though the "external" lag relations are 
perhaps the most significant determinants of rotating stall. The 
assumptions were: 
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(a) Straight entrance duct 

(b) Potential flow in entrance 

(c) Approximate form of relation between velocity components 
at entrance 

(d) Either straight or "suddenly infinite" exit duct 

(e) Inviscid flow in exit 

(f) Small-disturbance relation between velocity components in 
straight exit duct 

(g) Constant exit pressure in suddenly infinite exit duct 

The foregoing assumptions should be replaced by realistic relations 
obtained by experiment, with help from theory, treating realistic 
duct shapes. Velocity and direction, as well as hub and shroud static 
pressure should be measured close to the front and rear faces, but not 
too close in terms of individual blade effects. Extremely fast 
response will be needed to describe stall-zone boundaries which are 
quite abrupt (but not necessarily square waves) . 

Such measurements and theories will not only help rotating-stall 
prediction, but will show the way toward preventive measures, because 
it is in the entrance and exit flows that suppression of circumferential 
velocity can best be accomplished. 


CFD for Rotating Stall 


The comments in the previous paragraphs about the need to measure 
velocity relationships at the faces of the compressor imply the relevance 
of CFD. All the assumptions (a) -(g) require new calculations, certainly 
numerical, to apply theory to the overall compression system. Entrance 
"goosenecks" are not readily amenable to analysis, nor are the rather 
short diffusers of current practice. 


The Euler equations could be solved in these well-defined regions, 
subject to various disturbance wave forms for axial velocity. Perhaps, 
in addition to direct applications, it would prove possible to derive 
from the numerical results various kernel functions to generalize the 
velocity relations. 


Ultimately, we will be concerned with combined surge and rotating- 
stall theories. The two-dimensional, third order partial differential 
equations will surely require CFD treatment ultimately. Probably, the 
unsteady analysis of flow in blade passages is of less significance 
than that of entrance and exit regions. 
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APPENDIX F 


SYMBOLS 



C 

C 

P 

(C 

c 

D 

ID 

IE 

F 

F 

f 

G(g) 

g 

H 

h 


amplitude coefficient of disturbance (Eq. (2.6a) or (3.8)) 
stream-tube stress at entrance to IGV (Fig. 2.3) 
dimensionless reciprocal time lag, D/ (2 UNt) 

Fourier coefficient of g(0), Eq. (1.18a) 
denotes "break" of characteristic (Fig. 2.1) 
width of characteristic (Fig. 3.4(a)) 

denotes minimum of deep-stall characteristic (Fig. 2.1) 

disturbance of static-pressure rise coefficient across the 
compressor (Eq. (2,18)) 

rise of deep-stall line (Fig. 3.4(a)) 

constant, App. D 

mean wheel diameter 

slant of unstalled leg (Fig. 3.4(a)) 
slant of reverse-flow leg (Fig. 3.4(a)) 

steady pressure rise coefficient in blade passage, Eq. (1.3) 
height of drop into stall (Fig. 3.4(a)) 

stall propagation-speed coefficient, Eq. (1.11) and Fig. 1.4 

characteristic function, rotating stall absent (Figs. 3.1 and 
3.4(a)) 

axial flow coefficient perturbation in compressor, Eq. (1.4) 
height of entrance, App. D 

transverse velocity coefficient perturbation at entrance, 

Eq. (1.15) 

pressure-rise parameter at IGV entrance, Eq. (1.8) 
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parameter relating actual compressor-row lag to inertia 
lag of blade passage only 

depth of row in axial direction. Fig. 1.2(b) 

denotes maximum of characteristic at stall inception (Fig. 2.1) 
parameter defining lag tendency outside compressor, Eq. (1.34) 
number of stages in compressor core 

wave number, corresponding to number of stall cells 
pressure coefficient, Eq. (1.24) 
static pressure 

distance upstream from apex of wedge, App. D 

slope of characteristic at "break point" B, Fig. 2.1 

wheel speed at mean wheel diameter 

transverse velocity in laboratory frame 

average axial velocity 

axial velocity 

circumferential distance, in stall-zone frame, Eq. (1.5) 

distance normal to compressor face, in flow direction. 

Figs. 1.4,5 

new variable replacing h (Eq. (A8)) 

entrance angle, App. D 

throttle setting (Fig. 3.4(a)) 

stagger angle of blades. Fig. 1.2 

pressure recovery at entrance to IGV (Eq. (2.11)) 

difference of flow coefficient at points M and B of Fig. 2.1 

performance increase due to rotating stall (Eq. (2.5) or 
Fig. 3.1) 



E amplitude parameter of disturbance, Eq. (1.19) 

5 variable, Eq. (D4) 

rj nonlinearity parameter of disturbance (Eq. (2.6a)) 

9 angular location, in stall-zone frame, Eq. (1.5) 

X collection of terms analogous to oscillator mass, Eq. (1.39) 

p density 

a normal blade passage width. Fig. 1.2(b) 

T coefficient of pressure-rise hysteresis, Eq. (1.3) 

T* reduced time lag, Eq. (1.45) 

$ average flow coefficient, V/U 

(j) flow coefficient, v/U 

velocity potential, Eq. (1.13) 

Y upstream total to downstream static pressure-rise coefficient 

existing during operation 

compressor characteristic in absence of rotating stall. Fig. 1.1 
stream function, Eq. (1.23) 



caret: transformation defined by Eq. (3.11) 

double caret; after change of scale so that width of 
characteristic is 1 
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Subscripts 

cr critical or recovery value 

d design conditions 

e at compressor exit 

IV inlet guide vane 

OV outlet guide vane 

m amplitude of axial velocity disturbance 

ram conditions when axial velocity disturbance amplitude has 

its largest permitted value 

s far downstream static conditions 

T far upstream total conditions 

V combined inlet and outlet guide vanes 

0 just ahead of compressor entrance 

oo far upstream static conditions 
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